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1.0  SUMMARY 

1 .  1  Objective 

The  objective  of  this  contract  is  the  investigation  of  Saturable  Absorbers  at  10.  64 
and  3.  5m,  and  their  application  to  BALAD  receivers. 

1 . 2  Background 

In  the  BALAD  (Bleachable  Absorber  Laser  Amplifier  and  Detector)  concept,  the 
gaseous  absorber  is  designed  to  saturate  at  an  average  intensity  above  the  amplified 
spontaneous  emission.  Thus,  the  spontaneous  emission  and  other  background  light  is 
absorbed.  The  bleachable  gas  acts  au  an  optical  squelch  filter  suppressing  the  back¬ 
ground  radiation  which  does  not  coincide  with  the  signal  beam  in  direction,  polarization, 
frequency  and  time.  In  addition,  the  signal  beam,  being  of  somewhat  higher  intensity, 
is  focused  down  on  the  absorber,  "burns  a  hole"  in  it,  and  passes  through  with  relative¬ 
ly  little  attenuation.  Only  light  which  is  spacially  coherent  with  the  signal  gets  througl 
to  the  detector. 

During  the  first  year's  effort  under  Contract  F30602-7 1 -C -0024,  a  search  was 
made  for  useful  absorbers  in  the  10.  6m  and  3.  5m  ranges,  and  a  design  study  was  pur¬ 
sued.  As  a  result,  tits  gas  SF^  was  found  to  possess  adequate  properties  for  the  10.6|i 
wavelength  and  xenon  was  found  to  have  very  good  properties  at  3.  5m.  Demonstration 
of  the  BALAD  principle  at  3. 5 with  xenon  appeared  to  be  much  mere  straightforward 
and  considerably  less  costly  than  a  similar  experiment  at  10.6m  with  SF^.  The  present 
effort  is  a  logical  continuation  of  the  work  previously  performed  under  Contract 
F30602-7  1  -  C-0024,  and  includes  the  areas  defined  in  Tasks  4.  1.1  through  4.  1.3. 

1.3  Plan 

The  study  and  investigation  for  the  development,  design,  test,  and  analysis  of 
the  concept  of  a  Bleachable  Absorber  Laser  Amplifier  and  Detector  (BALAD)  has  been 
delineated  into  the  following  four  tasks: 

Task  4.  1.  1:  Development  of  design  techniques  for  the  BALAD  concept.  In  part 
this  will  involve  refinement  of  the  concepts  developed  under  Contract  F30602- 
7 1  -C-0024. 

Task  4.  1.2:  Develop  techniques  for  filtering  out  radiation  from  all  but  the  work¬ 
ing  line  in  the  10.6m  CO^  laser  amplifier  and  in  amplifiers  in  the  3.5m  spectral 
range  with  the  objective  of  determining  the  optimum  technique  to  be  jed  in  a 
BALAD  receiver. 

Task  4.  1.3:  Experimental  demonstration  of  a  iJALAD  receiver  utilizing  xenon 
both  as  the  absorbing  and  amplifying  gases;  in  particular,  using  a  pulse  modula¬ 
ted  xenon  oscillator  as  a  source,  measure  and  compare  with  theoretical  predic¬ 
tions  the  performance  of  the  BALAD  receiver  aj  a  function  of  input  pulse  energy, 
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duration  and  direction  of  propagation  with  reaped  to  the  parameters:  pulse 
energy  amplification,  amplification  bandwidth,  absorber  filter  bandwidth,  field 
of  view,  absorber  geometry,  no  detection  and  false  detection  error  rates  in  the 
pulsed  case,  signal  to  noise  ratio  in  the  CW  case. 

Task  4.  1.4:  System  Safety  Analysis.  In  the  conduct  of  design  and  fabrication 
of  hardware,  a  system  safety  analysis  shall  be  performed  and  documented. 

These  studies  shall  be  in  accordance  with  Paragraph  5.  8.  2.  1  of  MIL-STD-822 
where  applicable.’ 

1.4  Technical  Problems,  Methodology,  and  Results 

During  this  period  a  xenon-xenon  BALAD  receiver  system  has  been  operated 
and  tested.  The  results  are  reported  here.  The  results  of  the  theoretical  analysis  of 
the  propagation  of  a  realistic  optical  pulse  with  limited  transverse  extent  are  also  re¬ 
ported. 

The  theory  predicts  focut  sing  of  the  optical  pulse  in  a  saturable  absorber.  An 
analysis  is  presented  which  illustrates  how  this  effect  can  be  utilized  to  mprove  the 
BALAD  receiver  system.  The  following  subsections  summarize  the  developments  in 
each  of  these  topics. 

1.4.1  Experimental  Study  of  the  Xenon -Xenon  BALAD  Receiver  System 

A  BALAD  receiver  is  composed  of  a  laser  preamplifier  and  a  saturable  absorber. 
In  the  present  case,  both  components  employ  xenon  as  the  active  medium. 

The  desired  conditions  of  high  gain  with  a  largs  saturation  intensity  in  the  am¬ 
plifier,  compared  to  high  absorption  and  small  saturation  intensity  in  the  absorber, 
were  achieved  by  optimizing  the  operating  conditions  in  each  case.  In  the  1cm  I.D. 
amplifier  a  gain  coefficient  of  31dB/m  was  achieved  with  a  gas  mixture  of  2  Torr  helium 
plus  0.  1  Torr  xenon  and  a  current  of  20mA.  In  the  0.  5cm  I.D.  absorption  cell,  an 
absorption  coefficient  of  6CJB/m  was  achieved  with  low  satuiation  intensity  using  pure 
xenon  at  0.  1  Torr  and  a  current  of  125mA. 

The  experimental  investigation  of  the  performance  of  the  amplifier  and  absorber, 
both  separately  and  combined  together  as  a  BALAD  receiver  system  is  described  in 
Section  2.  4.  The  general  performance  was  in  agreement  with  the  theoretical  predictions 
However,  the  detailed  investigation  of  the  signal-to-noise  properties  of  this  BALAD 
receiver  system  waa  hampered  by  the  noise  generated  by  the  xenon  discharge  in  the 
absorption  cell.  We  were  unable  to  achieve  quiet  operation  of  this  discharge  under  the 
conditions  which  produced  the  strong  absorption. 


1>4*2  Transverse  and  Spatial  Evolution  of  Optical  Pulaea  in  a  Saturate 
Absorber  "  —  ’  - 

The  theoretical  analysis  of  the  propagation  of  an  optical  pulse  in  a  saturated 
absorber  has  been  significantly  advanced.  The  analysis,  which  is  described  in  Section 
3,  treats  the  realistic  case  of  a  pulse  of  limited  transverse  dimensions.  The  results 
indicate  that  a  strong  self  focussing  effect  occurs.  The  implications  of  this  effect  are 
summarized  in  the  next  subsection  and  are  discussed  in  Section  4.  The  propagation  of 
a  C  W  beam,  which  does  not  exhibit  the  self  focussing  effect,  is  discussed'in  Appendix  A 

1<4*3  Po^ible  Improved  Design  of  the  BALAD  Receiver  usine  Non-Linear 
Self-focussing 

The  existence  of  non-linear  self-focus  sing  suggests  the  possibility  of  using  geo¬ 
metrical  filtering  to  enhance  the  capabilities  of  the  BALAD  receiver.  In  the  BALAD 
receiver  as  presently  envisioned,  the  bleachable  absorber  must  attenuate  the  amplified 
background  noise  by  a  large  factor. 

As  detailed  in  the  previous  contract  reportc.1  this  required  the  use  of  a  rather 
long  absorber.  However,  as  detailed  in  Section  4,  non-linear  self-focussing  would 
allow  the  use  of  geometrical  filtering  to  accomplish  the  required  attenuatio...  The  sat¬ 
urable  absorber  would  the.,  function  primarily  as  a  wave  r,uide  although  the  attenuation 
which  it  proviaes  would  still  be  beneficial. 

1. 5  Implications  for  the  Department  of  Ocfe.:*e 

The  previous  study1  indicated  that  a  BALAD  receiver  for  10.6um  CC>2  laser 
signals  as  reported  in  RADC-TR-72-313,  is  feasible  for  application  in  an  optical 
radar.  Its  usefulness  would  be  to  detect  a  target  of  uncertain  position.  For  example, 
the  receiver  could  be  remote  from  the  transmitter.  The  BALAD  receiver  could  also 
sensitively  receive  several  separated  signals  simultaneously,  or  even  an  entire  image. 
The  spatial  and  spectral  resolution  of  a  SFfe  saturable  absorber  in  a  BALAD  receiver 
has  been  experimentally  verified. 

A  BALAD  receiver  utilizing  xenon  gas  could  function  in  a  shorter  range  but 
higher  resolution  optical  radar  than  the  CC>2  system.  This  wide-angle  receiver  also 
could  be  utilized  in  a  lightweight  secure  communication  system.  Discharge  noise  in  the 
xenon  absorber  may  limit  the  usefulness  of  th.»  xenon-xenon  BALAD  system. 

1  •  ^  Implications  for  Further  Research  and  Development 

The  difficulties  encountered  with  the  xenon  absorber  discharge  noise  point  up 
the  need  for  further  research  to  achieve  quiet  discharge  conditions  while  achieving  the 
desired  high  absorption  coefficient. 

The  self  focussing  effect,  indicated  by  the  numerical  •  nalysis  of  the  propagation 
of  a  spatially  limited  pulse  in  a  saturated  absorber  should  be  investigated  further.  If 
the  effect  -s  sufficiently  strong,  it  could  result  in  a  significant  improvement  in  BALAD 
receiver  design. 
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2.0  EXPERIMENTAL  STUDY  OF  A  XENON -XENON  BALAD  RECOVER  SYSTEM 
2.  1  Introduction 

We  have  studied  the  properties  of  the  BALAD  receiver  system  composed  of  an 
He-Xe  amplifier  and  an  Xe  absorber.  Since  the  pure  xenon  discharge  was  observed  to 
be  noisy,  in  most  of  the  measurements,  phase  sensitive  detection  was  used  to  improve 
the  signal-to-noise  ratio. 

2.  2  Equipment 

The  supportive  equipment  used  in  these  experiments,  grouped  by  category,  will 
be  briefly  described;  after  which,  a  more  detailed  discussion  of  the  elements  of  the 
BALAD  system  will  be  given. 

2.  2.  )  Power  Supplies 

Three  high  voltage,  low  ripple  power  supplies  capable  of  operating  in  the  range 
of  (i)  5kV,  100mA,  (ii)  20kV ,  50mA,  and  (iii)  lkV,  500mA  were  used  for  the  He-Xe 
l»ser,  He-Xe  amplifier,  and  the  Xe  absorber,  respectively.  A  Kepco  dc  supply  with 
range  10  V  and  6  A  was  used  for  the  hot  cathode  filament  of  the  Xe  absorption  cell.  Also 
a  Lambda  regulated  500  volt  dc  supply  was  used  to  bias  the  PZT  transducer  of  the  laser. 

2.2.2  Detectors 

A  liquid  nitrogen  cooled  (78  K)  Ge-Au  detector  was  used  to  detect  the  signal, 
ihc  bandwidth  of  the  detector  is  about  10  MHz.  The  output  from  the  detector  was  con¬ 
nected  to  a  Princeton  Applied  Research  JB5  lock-in  amplifier.  The  reference  frequency 
was  set  at  500Hz  which  was  the  frequency  of  the  chopper.  The  response  of  the  detection 
system  was  calibrated  against  an  Epply  thermopile.  The  laser  output  was  monitored 
with  a  Philco  L4530  InAs  detector  and  an  oscilloscope.  A  calibrated  Epply  thermopile 
with  sensitivity  6mW/mV  was  used  to  calibrate  all  detectors.  The  thermopile  voltage 
was  measured  with  a  Keithley  microvolt-ammeter. 

2.  2.  3  Oscilloscope 

In  all  the  measurements,  a  Tektronix  Type  551  dual  beam  oscilloscope  was  used 
ro  monitor  the  laser  signal  and  observe  the  noise  properties  of  the  system. 

2.  2.  4  Gas  Heating  and  Vacuum  Station 

The  station  consisted  of  a  diffusion  pump  and  a  mechanical  forepump.  One  liter 
bottles  of  helium  and  natural  xenon  from  Matheson  Gas  Products  were  attached  to  the 
Nation  oy  Granvillt -Pliilips  variable  leak  valves.  Gas  pressures  were  measured  with 
an  MKS  Baratron  Capacitance  Manometer.  The  plasma  tubes  of  .he  He-Xe  laser,  the 
He-Xe  amplifier,  and  the  Xe  absorber  were  attached  to  the  station  via  Granville -Philips 
ype  C  valves.  The  vacuum  station  could  be  pumped  down  to  10'7  Torr  after  thorough 
baking  of  the  glass  walls  and  the  valves  by  heating  tapes. 


2.2.5  Optical  Bench 

All  of  the  experiments  were  performed  on  a  98"x50''  heavy  cast  iron  bench. 

I  he  bench  was  isolated  from  the  floor  by  shock  absorbing  rubber  pads. 

2.3  Structure  of  the  Plasma  Tubes 

The  basic  elements  of  the  BALAD  system  were  the  xenon  amplifier  and  xenon 
absorber.  In  addition,  the  signal  source  was  a  xenon  laser. 

2.3.1  He-Xe  Laser 

The  He-Xe  laser  consisted  of  a  5  mm  I.D.  quartz  tube  with  a  dc  discharge  length 
of  26cm.  The  plasma  tube  had  2mm  thick  quartz  windows  oriented  at  Brewster's 
angle  (56  ).  The  quartz  windows  were  sealed  with  B.  F.  Goodrich  A-1177-B  epoxy. 

In  order  to  reduce  cathophoresis  of  xenon  during  the  dc  discharge,  a  long  non-exrited 
by-pass  was  provided  between  the  anode  and  the  cathode.  The  anode  was  made  of  nickel 
and  the  cold  cathode  was  made  of  aluminum  alloy  2024.  The  plasma  tube  was  mounted 
within  the  laser  resonator  which  had  three  invar  rods  held  in  circular  aluminum  bulk¬ 
heads.  The  resonator  mirror  mounts  had  three  point  fine  screw  orientation  adjustments. 
The  resonator  was  a  high-loss  Fabry-Perot  type.2  It  consisted  of  a  gold-coated  flat 
mirror  on  one  end  and  a  1mm  diameter  gold  dot  on  the  center  of  a  flat  quartz  substrate 
on  the  other  end.  The  gold  mirror  was  mounted  on  a  stack  of  PZT  elements  which  had 
a  response  of  53  A°/volt.  The  useful  laser  output  was  the  "diffraction  los  t »  around  the 
dot.  The  length  of  the  resonator  was  33.6  cm  long  which  corresponds  to  a  free  spectral 
range  of 

Av  =  ^-  =  447  MHz  . 

Due  *j  the  large  mode  pulling  effect,  the  actual  oscillation  range  of  the  laser  was  de¬ 
termined  to  be  only  239  MHz.  Since  the  full  doppler  linewidth  of  xenon  is  about  110MHz, 
single  mode  oscillation  was  obtained.  During  the  experiments,  the  He-Xe  la  er  was 
filled  with  a  pressure  of  3  Torr  of  helium  and  0.  08  Torr  of  xenon.  ’Vith  a  200  ’dl 
ballast  resistor,  quiet  operation  with  negligible  excess  noise  was  obtained  at  10rr.A 
discharge  current.  The  supply  voltage  was  3kV.  The  maximum  CW  output  from  t.ie 
laser  was  about  0.2mW. 

2.3.2  He-Xe  Amplifier 

Similar  to  the  plasma  tube  structure  of  the  He-Xe  laser,  the  amplifier  tube  con¬ 
sisted  of  a  1cm  I.D.  quartz  tube  and  a  dc  discharge  length  of  69cm.  Quartz 
Brewster  windows  were  used.  There  was  a  cold  aluminum  cathode,  a  large  ballast 
volume,  and  an  anode  to  cathode  bypass.  The  amplifier  was  filled  with  2  Torr  of 
helium  and  lOOmTorr  of  xenon.  At  a  current  of  20mA,  the  small  signal  gain  coefficient 
of  the  amplifier  was  31dB/m.  The  amplifier  plasma  tube  was  also  mounted  on  a  laser 
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structure  as  described  in  Section  2.  3.  1.  During  the  hole  burning  effect  observation  in 
xenon  absorber  as  will  be  described  in  Section  2.4.3,  the  He-Xe  amplifier  was  con¬ 
verted  to  a  laser  by  mounting  two  mirrors  to  form  a  high  loss  Fabry  and  Perot  reso¬ 
nator.  The  dot  reflector  had  a  1. 5mm  diamet  r  gold  dot.  The  maximum  output  power 
from  the  laser  wat  about  0.  6mW . 

2.3.3  Xe  Absorber 

The  absorption  cell  had  a  5mm  l.D.  and  a  dc  discharge  length  of  20  cm.  A  hot 
cathode  was  used  to  avoid  sputtering  at  high  currents  and  low  xenon  pressures.  The 
filament  was  made  from  a  12cm  long  piece  of  RCA  N91E1  oxide  coated  nickel  filament 
wire.  The  long  spiral  filament  was  supported  on  a  quartz  tube  to  avoid  sagging. 

Getters  were  included  to  aid  in  removing  impurities.  A  by-pass  was  also  provided 
between  the  anode  and  the  cathode.  Quartz  windows  were  used.  The  windows  were  only 
slightly  tilted  from  normal  to  the  optic  axis,  to  avoid  feedback  into  the  laser  without 
causi-.g  a  polarization  dependent  transmission  through  the  cell.  The  structure  of  the 
absoi  jtion  cell  is  shown  in  Fig.  1.  The  following  procedure  was  used  to  process  the 
hot  cathode:  The  entire  cell  was  first  outgassed  with  heater  tapes  and  pumped  down  to 
about  10‘6  Torr.  Then  the  filament  was  heated  by  enough  current  to  increase  the 
system  pressure  to  no  more  than  10-5  Torr.  The  heater  current  was  slowly  increased 
until  the  filament  was  heated  to  a  dull  orange  color.  The  getters  were  then  fired.  The 
whole  tube  was  baked  overnight  again  with  heating  tapes.  The  absorption  cell  was  then 
filled  with  xenon  as  desired  for  the  investigations. 

2.  4  Experimental  Study  of  the  Components  of  the  Xenon-Xenon  BALAD  System 


The  absorption  cell  was  first  systematically  studied  as  a  separate  system  ele¬ 
ment.  Th.  amplifier  was  similarly  studied  and  then  the  properties  of  the  combined 
system  were  investigated. 

2.4.1  The  Xenon  Absorption  Coefficient  Versus  Discharge  Current 

The  3.5  pm  signal  from  the  He-Xe  laser  was  chopped  at  500Hz  in  order  to  use 
phase  sensitive  detection  of  weak  signals  and  to  distinguish  fluorescence  from  the 
absorber  discharge.  A  50%  beam  splitter  was  used  to  direct  half  of  the  laser  power 
to  an  InAs  detector.  The  other  half  of  the  laser  output  was  sent,  single  pass,  through 
the  absorption  cell.  The  InAs  detector  signal  was  amplified  and  displayed  on  an 
oscilloscope  to  monitor  the  input  signal  power  and  control  the  oscillation  frequency 
near  line  center.  The  oscillation  frequency  was  manually  re-adjusted  to  the  «.toir' 
line  center  by  applying  dc  voltage  to  the  PZT  transducer.  The  output  signal  from  the 
absorption  cell  was  detected  by  a  liquid  nitrogen  cooled  Ge-Au  detector.  Tne  output 
from  the  detector  was  fed  into  a  lock-in  amplifier.  The  reference  signal  was  i  rived 
from  the  chopper.  A  plot  of  the  absorption  coefficient  (dB/m)  versus  discharge  cu  rent 


of  the  xenon  absorber  is  shown  in  Fig.  2.  It  can  be  seen  from  the  figure  that  maximum 
absorption  of  ~60dB/m  can  be  achieved  at  a  pressure  of  lOOmTorr  and  at  a  current  of 
125mA.  Comparing  this  to  previous  data3  one  finds  that,  quite  analogous  to  the 
gain  coefficient,  the  maximum  absorption  coefficient  is  inversely  proportional  to  the 
tube  diameter. 

2.4.2  The  Absorption  Coefficient  Versus  Input  Power 

Our  analysis  of  the  saturation  of  absorption  coefficient  vs.  input  power  is  based 
on  the  work  of  C.P.  Christensen  et  al.  f  who  studied  the  gain  versus  input  power 
characteristic  in  a  C02  amplifier.  The  absorption  coefficient  for  an  inhomogeneously 
broadened  medium  is  given  by 


where 

oq  is  the  small  signal  absorption  coefficient 
I  is  the  intensity  at  z 

Ig  is  the  saturation  intensity  of  the  medium. 

For  small  intensity,  i.e.,  l/lfl  «  1,  Eq.  (1)  can  be  approximated  by 

1  dl  u  1  I  %  I;) 
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After  integrating  this  equation  over  the  length  of  the  absorber  L,  the  result  is 

t  »  -o  L  . 
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out  s 

Equation  (3)  serves  as  a  means  to  determine  Is  in  a  plot  of  ^n/^out  VB- 

If  the  laser  beam  has  a  constant  size  and  shape  throughout  the  absorber,  the 

intensity  can  be  integrated  over  the  beam  cross-section  to  yield  a  relation  between  the 

input  and  output  power.  For  example,  if  the  incident  beam  iB  Gaussian  of  the  form 

I.  =  I  exp(-r  /w  )  where  w  is  the  beam  radius  at  l/e  intensity,  then  the  power  input 
in  o  r  ' 

is  given  by  integrating  over  the  area  of  the  beam,  i.e.  , 
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Absorption  coefficient  of  5mm  I  D.  Xenon 


Similarly,  the  output  intensity  is  given  by 
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Substituting  these  results  into  Eq.  (3)  and  integrating  over  the  cross-sectional  area 
yields : 
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Equation  (6)  shows  a  straight  line  in  the  P  /P  vs-  P  plot.  The  intercept  at 

(P  /P  )  =  0  gives  the  value  of  saturation  intensity  I  which  is 
in'  out  **  s 


iru 


(7) 


A  typical  plot  of  the  ratio  of  input  power  to  output  power  vs.  input  power  is  hown 
in  Figs.  3a  and  3b.  From  Fig.  3b,  the  slope  at  the  small  input  power  gives  an  inter¬ 
cept  at  l6gt  W .  The  average  beam  size  in  the  absorber  was  about  u>=  0.  5mm  and 
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From  Expression  (7),  the  saturation  intensity  of  the  absorber  is  estimated  to  be 

1  *  0.  51mV.f/cm^  in  good  agreement  with  the  previous  measurements  of  xenon  gain  sat- 

uration. 

2.4  3  The  Output  Power  Versus  Frequency 


The  3.5pm  output  from  a  strong  He-Xe  laser  was  sent  through  the  absorption 
cell  after  being  collimated  by  a  pair  of  concave  mirrors  with  radii  of  curvature 
Rj  =  92cm  and  =  42cm  used  as  a  si  nple  telescope.  A  weak  signal  from  a  set  ond 
He-Xe  laser  which  served  as  the  probe  signal  was  also  directed  through  the  above  col¬ 
limator  before  reaching  the  absorption  cell.  The  weak  beam  was  arranged  to  overlap  the 
strong  beam  within  the  absorption  cell,  but  it  propagated  at  a  small  angle  relative  to 
the  strong  beam  so  that  the  beam*  could  be  separated  and  only  the  weak  signal  rc;u  bed 
the  Ge-Au  detector.  The  frequency  of  the  strong  laser  output  was  fixed.  The  P7  1 
transducer  of  the  weak  laser  was  sw  *pt  by  a  sawtooth  output  voltage  .rom  an  oscilloscope 
to  obtain  the  frequency  tuning  profile  of  its  output  power.  A  large  amount  of  noise  from 
the  absorption  cell  was  obtained  along  with  the  weak  signal  when  the  cell  was  operating 
at  a  current  of  100mA  and  at  a  pressure  of  lOOmTorr.  It  was  impossible  to  observe 
the  hole  burning  effect  directly  from  the  oscilloscope.  Phase  sensitive*  detection  was 
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used  10  improve  the  signal-to-noise  ratio.  The  weak  signal  wac  chopped  at  500Kz  and 
sent  through  the  absorption  cell.  The  detected  signal  then  went  to  a  iock-in  amplifier. 
The  output  from  the  lock-in  amplifier  was  connected  to  tiie  oscilloscope  which  was 
swept  synchronously  with  the  weak  laser.  The  hole  burning  effect  was  observed  as  show 
in  Fig.  4. 
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The  absorption  coefficient  of  the  weak  signal  in  the  presence  of  a  strong  satura¬ 
ting  beam  of  intensity  I  at  frequency  is  given  by 

oo 

a[v,I(v1)]=  J  o(v^)AN[v^,  Ifv^Jdv^  , 
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where 


0(v^  )  = 
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AN  v^,  I(Vj)] 


is  the  population  difference  between  the  lower  and  upper  state  of 
the  atom  in  the  presence  of  a  strong  saturating  signal, 


v 


o 

Av 


is  the  atonrc  line  center  of  t.nse  atoms  with  non-zero  velocity, 

is  the  atomic  line  center  of  those  atoms  with  zero  velocity, 

is  the  full  homogeneous  line  width  of  the  atom  at  half  intensity, 

■r  =  {>./ 2tr)  (l/t^v)  is  the  absorption  cross-section  of  the  atom, 

\  is  the  wavelength  of  the  laser, 

is  the  spontaneous  lifetime  of  the  atom. 


11  we  neglect  the  fluctuation  of  the  population  due  to  the  beat  frequency  between 
the  tw'o  lasers,  the  population  difference  can  be  written  as 
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where 


AN 
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is  the  population  difference  of  the  medium  without  incident  field, 
is  the  saturation  intensity  of  medium, 


P(v^) 


is  the  weighting  funct  on  of  the  atom  at  with  respect  to  vq. 


Hie  1  ittcr  quantity  is  given  by 
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Weak  muni  Iransmi-'  ion  ol  the  xenon  absorber  while 
!'  is  .alnrated  by  .1  strong  sien.il  \s.  short  iavity  Uini>.u 
length.  Horizontal  s>nie  i~  .a;>j>roximately  t*.  K-l  MHz/ box 


where 
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=  / 
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is  the  full  doppler  line  width, 


k  =  Boltzmann's  constant,  T  =  temperature,  M  =  atomic  weight,  and 

c  =  velocity  of  light. 

Substituting  Eqs.  (9-11)  into  Eq.  (8)  and  performing  the  integration  in  the  limit  of  broad 
doppler  line  width,  one  finds  , 
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is  the  small  signal  absorption  coefficient  of  the  medium.  From  Eq.  (12),  we  have  at 
V  =  V, 
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The  output  intensity  IQut  of  the  weak  signal  from  an  absorber  of  length  L  is  related  to 
the  input  intensity  1.^  by 

IoutsIincxP<’oL)  <15> 

v'here  Q  is  given  by  (12).  At  v  =  Vj,  the  weak  signal  output  from  the  absorber  in  the 
presence  of  the  strong  signal  can  be  obtained  from  (15),  i.e.,  1^  =  1,^  exp!"-a(  Vj)L]  r 
Iin  exP  ‘Q0( Vj)L/^^ 1  +  (l/ls)  ].  When  there  is  no  strong  signal,  the  weak  signal  output 
from  the  absorber  at  v=  is  simply  given  by  IQut  =  exp  Vj  )L].  The  depth  of  the 
hole  is  thus  given  by  the  difference 

Ii"Cexp(-^fe)'exp,-a»L,J  • 


The  full  width  of  the  hole  at  half  of  the  hole  depth  can  be  obtained  by  solving  the  follow  - 
ing  equation 
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which  gives  after  some  calculation 
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The  experimental  data  gives  at  v  =  , 
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Hence  one  finds  (i/i g )  «  0.33.  From  (16),  v.c  then  obtain 


2 (v  -  Vj)  *  l.  97  Av  .  (17) 

The  full  natural  line  width  of  xenon  at  3.  5tiin  is  about  3. 66  MHz,5  and  the  pres¬ 
sure  broadening  of  the  homogeneous  line  width  is  about  60MHz/Torr.3  For  lOOmTorr 
of  xenon  used  in  the  absorber,  the  pressure  broadened  homogeneous  line  width  is 
-6  MHz.  Hence  we  have  Av  9.  75  MHz.  From  Eq.  (17),  we  can  estimate  the  full 
width  of  the  hole  burned  in  the  absorption  curve  is  about  19  MHz  which  is  in  fair  agree¬ 
ment  with  the  measured  16  MHz  width. 

2.  5  Transmission  Characteristics  of  the  BALAD  System 

The  3.  5bm  output  from  the  He-Xe  laser  was  sent  through  the  He-Xe  amplifier 
after  reflection  from  a  long  radius  concave  mirror  of  R  =  500cm.  The  amplified  beam 
was  focused  into  the  center  of  the  Xe  absorber  by  a  short  radius  concave  mirror  of 
R  =  41.8cm.  The  output  from  the  absorber  was  collected  by  a  concave  mirror  of  R  =  30im 
onto  a  Ge-Au  detector  at  77°K.  Phase  sensitive  detection  was  used  in  the  measurements 
as  described  in  Section  2.4.  1.  A  beam  splitter  was  used  to  monitor  part  of  the  laser 
output.  The  schematic  diagram  of  the  experimental  set  up  is  given  in  Fig.  5.  Part  of 
the  signal  from  the  Ge-Au  detector  was  fed  to  the  oscilloscope  to  observe  the  signal  and 
noise  properties  of  the  BALAD  system.  Glass  microscope  slides  were  used  to  attenuate 
the  laser  output  power.  The  slides  have  an  attenuation  factor  of  about  50  percent. 
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Fi«{.  S  -  Experimental  -el-up  for  measurements  of  transmission  eh.  -aclensUc.  of  a  Xe-Xe  BALAD  system 


The  following  procedures  were  followed  during  the  measurements.  The  signal 
power  from  the  He-Xe  laser  was  measured  by  the  lock-in  amplifier  with  the  discharge 
in  both  the  amplifier  and  absorber  off.  The  noise  voltage  of  the  signal  was  mcajureo 
from  the  oscilloscope.  Then  the  amplifier  was  turned  on.  The  amplified  signal  power 
and  the  noise  were  measured  by  the  lock-in  amplifier  and  the  oscilloscope.  Finally, 
bo*h  the  amplifier  and  the  absorber  were  turned  on  and  the  transmitted  signal  power 
and  noise  were  measured.  The  glass  attenuators  were  inserted  in  front  of  the  laser  on? 
at  a  time  and  ths  above  procedure  was  repeated.  During  the  experiment,  the  current 
of  the  He-Xe  laser  was  fixed  and  the  laser  oscillation  frequency  was  readjusted  to  the 
line  center  by  manual  control  of  the  laser  PZT  transducer  voltage.  The  amplifier  pres¬ 
sure  was  2Torr  of  helium  and  lOOmTorr  of  xenon.  The  discharge  current  was  main¬ 
tained  at  20mA.  The  xenon  pressure  in  the  absorber  was  lOOmTorr  and  the  discharge 
current  was  maintained  at  about  125mA.  The  plot  of  transmission  (=  output  power/ in¬ 
put  power)  vs.  input  power  for  the  amplifier  and  the  BALAD  system  are  shown  in  Figs. 
6a,  6b,  and  7.  In  the  case  of  absorber,  the  ratio  input  power/output  power  vs.  input 
power  is  plotted  in  Figs.  3a  and  3b. 

For  the  high  helium  and  xenon  pressure  used  in  the  amplifier,  the  gain  becomes 

essentially  homogeneously  broadened.^  The  saturation  intensity  of  the  amplifier  can  be 

q 

roughly  determined  fcy  the  intercept  of  the  slope  at  small  input  power  by  the  relation 

iru> 


From  Fig.  6b,  one  finds  *  12.6pW, 
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and  the  beam  size  in  the  amplifier  is  about  u  =  3mm.  The  saturation  intensity  of  the 
amplifier  is  calculated  from  Eq.  { 1 H)  to  be  about  2.5mW/cm  . 

Figure  7  shows  a  peak  in  the  transmission  of  the  signal  through  the  BALAD 
system  which  occurs  at  input  power  to  the  amplifier  of  about  0.28p\V.  Since  the  satura¬ 
tion  intensity  of  xenon  absorber  is  very  small,  saturation  processes  started  when  the 
signal  power  to  the  amplifier  was  onl^  0.02mW.  At  such  an  input  power,  the  amplifier 
was  not  saturated.  The  input  signal  was  amplified  by  a  factor  of  one  hundred  and  thirty 
times  through  the  amplifier  to  give  ar  input  of  2.6i.  \V  to  the  absorber.  At  the  above 
input  |  'wer,  there  was  only  about  a  factor  of  14.  5  attenuation  through  the  absorber  is 
compared  to  the  factor  of  16  unsaturated  small  signal  ittenuation.  Hence  there  was  i 
relative  increase  of  signal  power  through  the  BALAD  system.  As  the  input  «ignil  level 
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to  the  amplifier  was  increased,  i*  was  linearly  amplified  through  the  amplifier  and  the 
absorber  was  nonlinearly  decreasing  its  absorption  ability  wliile  it  continued  to  I  c 
saturated.  More  transmission  through  the  BALAD  system  would  be  expected  When 
the  input  signal  to  the  amplifier  reached  the  level  of  about  O.ZHuW,  it  started  to  satu¬ 
rate  the  amplifier.  The  amplification  factor  decreased  while  the  absorber  was  almost 
at  a  constant  attenuation.  The  transmission  through  the  BALAD  system  dec  re  ised. 

As  the  input  signal  to  the  amplifier  was  further  increased,  the  amplifier  was  saturated 
faster  than  the  absorber  as  can  be  seen  from  Figs.  3a  and  ba.  The  transmission 
continued  to  decrease  is  shown  in  Fig.  7. 

In  order  to  investigate  the  propagation  of  the  beam  through  the  BALAD  system, 
one  has  to  start  from  the  equations  which  govern  the  propagation  of  Gaussian  beam. 
The  following  discussion  will  presume  "free  space"  propagation.  The  modifications 
due  to  radial  variations  in  gain  or  absorption  will  be  discussed  in  Appendix  A.  The 
evolution  of  short  pulses  will  be  discussed  in  Section  3  and  Appendix  B.  For  a  given 
waist  spot  size  uiq  at  z  =  0,  the  spot  size  w  and  the  radius  of  curvature  of  the  be  mi ,  R , 
a  distance  z  are  given  by 


w(z)  =  “qV  1  +  (>z/iru-o) 


R(z)  =  z  Id  +  (it uT/Xz) 


2, 


(19) 

(20) 


where  \  is  the  wavelength  of  the  beam.  Solving  Eqs.  (19)  and  (20)  for  and  z  in  terms 
of  c.(z)  and  R(z),  one  gets 


z  = 


R 


IM-^V 


(21) 


v 1  + 


Wher  a  Gaussian  beam  with  radius  of  curvature  Rj  incidents  on  a  concave  mirror  with 
local  lerg  h  f.  The  radius  of  curvature  R of  the  reflected  beam  will  be  transformed 
acrordLig  to  the  lens  law 


_L  _L  i 

R2_K1  ‘f 


(23) 


The  radius  of  curvature  of  the  beam  from  the  3.  5i-m  He-Xc  laser  .it  mirror  M  , 

L. 

was  about  2n5cm  and  the  spot  size  at  M,  was  found  to  be  -  3mm.  From  F.qs.  (21), 

u  2 

(22),  and  (23),  the  waist  spot  size  w  and  its  distance  Lj  from  M,  is  calculated  loin 
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about  2.8mm  and  290cm,  respectively.  The  angular  field  of  view  of  the  amplifier  for 
1  mode  is  given  by 


n 

am 


1  mode 

P 


)  = 


0.  5x  10'6 


ste  radians 


(2-1) 


The  angular  field  of  view  of  the  amplifier  in  all  modes  is  given  by 

^amp(a"  modes)  =  ^  =  0.  24x  10  ^  steradians  (29) 

where 

Dj  =  1  cm  is  the  diameter  of  the  amplifier 

dj  =  llhcm  is  the  distance  from  the  mirror  M?  to  the  window 
of  the  amolifier.  ^ 


The  number  of  modes  (or  resolution  elements)  m  of  the  amplifier  is  given  by 

n  (all  modes) 

m  (amplifier)  =  ^  mpd  =  48  .  (2b) 

amp  ' 


By  substituting  z  =  M£M3  +  M3M4  -  =  240  +  1  10  -  290  =  60cm  into  Eq.  (19)  and  Eq. 

(20),  one  finds  the  spot  size  and  radius  of  curvature  of  the  beam  at  mirror  to  be 
“M4  =  2.8  mm  and  R  =  8.6  xlO5  cm.  Using  Eqs.  (21),  (22),  and  (23),  the  waist  spot 
size  and  its  distance  from  M4  is  given  by  =  0.  17  mm,  and  L2  *  20.  9cm.  The 
angular  field  of  view  of  the  absorber  for  1  mode  is  given  by 


K  -  4 

^abs^*  111  °de)  =  - 7  =  1.35x10  steradians 


(27) 


iru; 


The  angular  field  of  view  of  the  absorber  in  all  modes  is  given  by 


tr(D2/2) 


^abs^3^  rood®8)  =  - - - j  =  1 . 97  x  10  ^  steradians 


(L2-d2) 


(28) 


where 


D2  =  0.  5cm  is  the  diameter  of  the  absorber 

d2  =  10cm  is  the  distance  from  M  to  the  window  of 
the  absorber. 


The  number  of  modes  m  of  the  absorber  is 


m  (absorber) 


nabs(all  modes) 
^abs(1  modc)~ 


(2‘1) 
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Ihe  spot  size  at  mirror  M$  is  about  ^  *  1.3cm  and  the  radius  of  curvature  of  the 

beam  after  reflecting  from  is  about  26cm.  Similarly,  from  Eqs.  (22)  and  (21), 

the  waist  spot  size  and  the  distance  is  calculated  to  be  0.  22  x  10*2cm  and  26  c  m, 

respectively.  The  number  of  modes  that  the  Ge-Au  detector  is  able  to  receive  is  Jive, 
by 

-  Aiea  of  detector  7.9x  10'3cm2  ... 

2  -  - 5 - 500  .  (3m 


2.  6  Conclusions  of  the  Experimental  Study 

This  completes  the  description  of  the  experimental  demonstration  of  the  xenon- 

xenon  BALAD  system.  The  conclusions  of  this  study  are  presented  along  with  other 
results  in  Section  6. 


3.0  THEORETICAL  INVESTIGATION  OF  COHERENCE  PULSE  PROPAGATION  AND 
SELF  FOCUSSING  OF  COHERENT  PULSES 

3.  1  Introduction 

When  optical  pulses  propagate  coherently  through  resonant  media,  dramatic  dif¬ 
ferences  in  behavior  from  ordinary  propagation  occur.  These  effects  depend  on  the 
resonant  interaction  of  the  polarizable  medium  with  pulses  short  compared  to  tlu*  phase 
memory  time.  An  example  is  the  phenomenon  of  "Sell-Induced  Transparency."  A 
review  paper  by  G .  L.  Lamb,  Jr .  ^  surveys  the  work  in  this  area .  Gibbs  and  Slushc  r;’  ** 
have  obtained  experimental  verification  of  important  aspects  of  coherent  pul  it  propa>  >- 
tion  phenomena.  In  the  published  work,  the  emphasis  has  been  on  configurations  with¬ 
out  spatial  variations  in  the  directions  transverse  to  the  direction  of  propagation. 

Coherent  propagation  with  transverse  variation  has  been  discussed  by  McCall  and  H.ihn^ 

using  a  generalization  of  one -dimensional  "Pulse  Area-Pulse  Energy"  theorems. 

1 2 

Barone  and  Chi  have  shown  how  the  inclusion  of  transverse  effects  in  the  theory  can 
lead  to  qualitative  changes  in  ’he  behavior  of  the  predicted  steady  state  pulses.  In  i  on  - 
nection  with  studies  invoicing  an  aosorber  driven  by  strong  pulses,  Rhodes  and  Szoke'  * 
have  noted  efffets  which  they  attributed  to  transverse  variations  in  the  optical  beam 

Transverse  phenomena  have,  of  course,  been  emphasized  in  the  work  dealing 
with  self-focussing  of  optical  pulses,  but  these  usually  have  been  concerned  with  the 

effects  of  non -res on.mt  non -liima rities  .  A  review  paper  by  S.A.  Akhmanov,  A.P. 

1  4 

Sukhorukov,  and  R.V.  Khokhlov  contains  a  bibliography  of  work  up  to  19b7.  Self- 
focussing  involving  resonant  interactions  has  been  studied  by  A.  Javan  and  P.  L.  Kelley  *  ^ 
and  by  D.  Grischkowsky  and  J.A.  Armstrong.^  Both  of  these  investigations  rel>  on 
the  dispersive  part  of  the  susceptibility  on  the  wings  of  the  resonance  line  to  provide 
the  non-linear  index  of  refraction.  The  response  of  the  medium  is  describable  in  terms 
of  a  simple  intensity  dependent  susceptibility,  rather  than  a  more  general  functional  of 
the  applied  field.  In  the  first  paper,  this  was  possible  because  the  field  is  slowly  vary¬ 
ing  in  a  relaxation  time,  and  in  the  second  because  its  variation  in  time  is  consistent 
with  the  approximations  which  lead  to  "adiabatic  following.  " 

When  there  is  a  susceptibility,  the  mathematical  formulation  of  self -foe  us  sing 
can  be  made  equivalent  to  the  mechanical  problem  of  particle  motion  in  a  potential 
well  with  friction;  the  shape  of  the  well  being  determined  by  the  form  of  the  susceptibility. 
Thus,  one  can  easily  visualize  the  formation  of  self-focussing  modes.  In  the  more 
general  case,  this  useful  model  is  not  available  and  greater  reliance  is  had  on  nutneri  ul 
integration  of  the  coupled  field  material  equations. 

J.A.  Fleck,  Jr.  and  C.  Layne*'  have  included  the  effect  of  resonant  interactions 
in  their  numerical  study  of  self-focussing  damage  in  a  Nb-glass  amplifier.  Fluir  for¬ 
mulation  was  capable  of  handling  coherent  effects  sim  e  the  resonant  interaction  \v  is 
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described  in  terms  of  Bloch's  equations.  However,  for  the  system  of  interest,  these 
equations  were  solved  in  the  rate  equation  approximation  and  hence  did  not  include 
coherent  effects. 


We  have  been  studying  the  propagation  of  short  optical  pulses  through  resonant 
media  with  emphasis  on  transverse  effects  associated  with  coherent  interactions  be¬ 
tween  the  field  and  the  material  system.  When  the  incident  field  has  substantial  trans¬ 
verse  variation  over  the  input  aperture  significant  alterations  of  transverse  structure 
are  to  be  expected  under  certain  conditions,  even  within  the  "near  field"  region  of  the 
effective  input  aperture.  This  expectation  is  based  on  the  dependence  of  the  velocity  of 
a  one-dimensional  pulse  on  the  pulse  area  and  shape.  The  dominant  feature  of  coher¬ 
ent  pulse  propagation  in  absorbers  are  the  reshaping  toward  a  steady  state  form  a  id 
pulse  speeds  that  are  considerably  less  than  the  speed  of  light  in  the  medium.  In  the 
three-dimensional  problem,  we  would  expect  a  region  surrounding  each  ray  to  at  first 
propagate  like  the  corresponding  one -dimensional  input  field  of  the  same  amplitude. 
However,  after  a  certain  distance  longitudii  al  reshaping  will  lead  different  rays  to 
carry  pulses  with  different  speeds.  This,  in  turn,  will  lead  to  transverse  reshaping. 
When  this  occurs,  the  radial  derivatives  in  the  wave  equation  can  lead  to  significant 
transverse  phase  variations  and  energy  flow.  Because  of  the  complexity  of  the  coupled 
field-material  system  equations,  we  have  used  numerical  integration  to  investigate 
these  effects. 


3.  2  Field  and  Material  System  Equations 

In  our  work,  the  electric  field  is  assumed  to  be  plane  polarized  transverse  to 
the  direction  of  propagation,  z,  and  of  the  form 

E=5c+fe  ,  (3 

with 


0  =  wt  -  unz/ c  ,  (32) 

where  n  is  the  index  of  refraction  associated  with  the  non-resonant  part  of  the  material 
system  and  w  is  the  carrier  frequency.  The  resonant  part  of  the  material  polarization 
is  taken  to  be  of  the  form 


P  = 


-if 


-iO 

e 


+ 


(33) 


Assuming  that  the  complex  amplitudes  if  and  f  vary  slowly  in  their  z  and  t  dependence 
compared  to  the  high  frequency  variations  of  the  carrier,  the  wave  equation  in  cylindri¬ 
cal  coordinates  simplifies  to 


L'1 


12.0,  no  „  2ttu,,_ 

Cep  cz  c  2  j  nc 

c. 


(3  1) 
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This  form  of  the  wave  equation  has  been  used  in  discussions  of  self-focussing1  *  with 
the  polarization  amplitude  G  expressed  as  a  non-linear  function  of  the  field  amplitude 
3  ■  For  coherent  interactions,  the  polarization  amplitude  is  a  functional  of  the  lield. 

In  order  to  be  capable  of  describing  coherent  effectB,  we  model  our  material  system  by 
a  gas  of  two  level  atoms.  A  general  description  of  such  a  system  would  include  the 
effects  of  degeneracy  and  inhomogeneous  broadening.  The  polarization  amplitude  would 
be  a  sum  over  i  degeneracy  index  j  and  an  integral  over  an  inhomogeneous  line 

P  =  ^  N  /  g(Ai*)  P  (Au)dAu  (3r>) 

J  J  J 

where  is  the  number  density  corresponding  to  the  generacy  index  j  and  g(Au,)ri(A«,) 
is  the  fraction  of  molecules  in  the  interval  dAu  whose  central  resonant  frequency  differs 
from  the  inhomogeneous  line  center  by  the  frequency  interval  Ac*.  Each  polarization 

g 

amplitude  P.(Acj)  satisfies  the  Bloch  equations.  \ 

J  \ 

d  i  Kl2 

(3t  +  T^+  lA“)pj  =  ~ 5  wj  ( 

(£  +  T")(Wi  -  W°)  8  -  ‘  +  /  P)  .  (37) 

1  1  J  J  <1 

In  this  section  we  will  report  on  the  study  of  transverse  effects  in  non-degenerate 
homogeneous  systems.  In  Eq.  (35),  the  degeneracy  index  j,  takes  on  only  the  varne  1, 
and  the  line  shape  function  g(Aw)  is  a  delta  function,  g(A<*)  =  g(Aw).  We  have  studied 
amplifiers  and  absorbers  and  cases  where  the  input  pulse  duration  is  sho^  t  ns  well  as 
comparable  to  the  relaxation  times.  When  the  pulse  duration  is  comparable  to  and 
longer  than  the  relaxation  time,  it  is  appropriate  to  describe  the  system  in  terms  ol  . 
linear  gain  or  loss  coefficient  a  in  units  of  cm  In  our  notation  a  is  given  by 


Q  =  NW°  iiii-  1  .  (jb) 

o  nc  ft  2  '  ’ 

This  quantity,  the  Fresnel  number  of  the  aperture  (diameter  2p  )  corresponding  to  the 
distance  oq  ,  PQAQ0  ,  and  ti  e  ratio  of  relaxation  times  T^/T^  provide  three  independ¬ 
ent  parameters  to  define  the  solution.  We  have  chosen  these  parameters  to  be  con 
si  stent  v/ith  X  =  10  cm,  qq  =  ±(l/20cm),  Pq  =  1cm,  =  10'9sec,  Tj/T^  2.  The 
sign  of  qq  depends  on  whether  we  are  dealing  with  an  attenuator  or  absorber. 

When  the  pulse  duration  T  ,  is  short  compared  to  the  relaxation  time  T-.,  it  is 

P  2’ 

no  longer  appropriate  to  refer  to  the  linear  absorption  length.  In  that  case,  we  scale 

to  Ip  the  quantities  that  were  previously  referred  to  T,,  and  take  T  =  10"  *  sec, 

T  /T  »  1.  In  the  calculations  reported  on  below,  we  use  T.,  =  )  0T  ,  and  I--,  1  . 


In  this  section,  we  v;ll  describe  the  results  of  our  calculations  of  transverse- 
effects  associated  with  coherent  propagation  through  absorbers.  We  will  relate  the 
results  of  our  numerical  analysis  of  propagation  in  j  dimensions  with  corresponding 
calculations  for  the  one-dimensional  configuration.  It  will  be  seer,  mat  significant 
self-focussing  can  occur  when  coherent  pulses  propagate  through  a  resonant  absorbing 
medium.  We  differentiate  between  "transverse  reshaping"  and  "self-focussiv  g .  "  The 
latter  is  associated  with  the  transverse  flow  of  energy.  This  "current"  may  be  identi¬ 
fied  from  the  energy  equation  derivable  from  Eqs.  ( 34) -( 37) . 


2  i  i  2 

7. 1.-4  n-^~ 


Nhtj  ...  Nhu  [ 

+  TwJ‘Tl 


...  .w-w 

Nhu  /  o 


where  the  radial  and  longitudinal  components  of  the  intensity  vector  ,  and  J z  are  given 


J  =  £-i(<JVf*  -  s\s) 

p  2u 


Jz  =  nc|<?| 

In  terms  of  the  amplitude  |  S  |  and  phase,  c? •  of  the  complex  amplitude: 


S  =  |  8  |  e1 


we  have 


JP- ir  1^1  ip  • 

Thus  transverse  energy  flow  is  associated  with  a  radial  variation  of  the  phase  of  the 
complex  amplitude  <S . 

3 . 3  Transverse  Effects  in  Absorbers 

3.3.1  One  Dimension 

In  our  discussion  of  absorbers,  we  first  consider  the  case  of  one -dimensional 
pulse  propagation  with  the  pulse  duration  much  shorter  than  both  Tj  and  T^.  Analytic- 
distortionless  solutions  of  the  material-field  equations  are  known8  of  the  form 

i-sech  H3I 

*  TP  TP 

where  v  is  the  pulse  velocity  given  in  'erms  of  the  pulse  length  by 


!=  a  +  _2_  T2 

v  c  T,  p 


2B 


This  is  the  2ir  pulse: 


-00  h 


(45) 


In  Fig.  8,  we  present  the  results  of  numerical  integration  of  the  one -dimensional  field 
equation  for  the  case  that  the  input  pulse  is  a  gaussian  in  time  of  input  area  2ir. 

Figures  8(a)  through  8(c)  show  the  pulse  reshaping.  In  this,  and  succeeding 
figures,  the  symbol  t  stands  for  the  time  counted  from  the  instant  when  a  signal  travel 
ling  at  the  speed  of  light  in  the  medium  reaches  the  point  z,  thus 


1  he  units  °i  T  are  nanoseconds.  After  300cm  the  pulse  shape  closely  matches  the  sech 
form.  From  the  measured  width  after  300cm  the  analytic  expression  for  the  pulse 
velocity  is  (v/c)  =  0.7.  In  Fig.  8(d),  we  have  plotted  the  integrated  intensity 
(Energy/cm  =  /d  T  |  <?  |  ),  pulse  area,  and  time  delay  of  the  peak  all  against  the  propa¬ 
gation  distance  z.  The  time  delay  is  the  difference  between  the  final  and  initial  values 
of  t  at  the  peak  of  the  pulse.  This  quantity  measures  the  velocity  of  the  peak  relative 
to  the  speed  of  light. 


c  Lc(A  r/Az)  +  1  J  ’  (47) 

where  At  is  the  time  delay  accumulated  in  the  interval  Az.  We  note  that  the  pulse  area 
remains  at  the  inpvt  value  2tt  and  the  pulse  energy  decreases  linearly  with  distance  as 
the  pulse  reshapes  from  gaussian  toward  sech.  This  decrease  is  associated  with  the 
non  -vanishing  (but  small)  value  of  T^.  The  pulse  peak  mo\es  more  slowly  than  the 
speed  of  light.  After  300cm,  the  numerical  pulse  velocity  is  0.7  times  the  speed  of 
light.  This  corresponds  to  the  speed  of  the  sech  pulse  of  the  same  width. 

3.3.2  Three  Dimensions 

Figure  8  is  to  be  contrasted  with  Fig.  9  which  presents  the  results  of  numerical 
integration  of  the  three-dimensional  field  equations  for  the  same  parameters.  The  in¬ 
put  pulse  has  a  gaussian  variation  with  radius,  p,  as  well  as  with  time,  T,  in  the  input 
plane  z  =  0.  This  is  illustrated  in  Fig.  9(a)  where  the  units  of  p  arc  cm.  In  this  figure, 
the  dependence  of  the  field  amplitude  |  f  |  on  t  is  shown  for  three  values  of  p  ,  p  =  0, 
Vl4cm,  and  s/7cm.  Up  to  about  40cm,  the  variation  with  t  along  the  axis  f  =  0  is 
the  same  as  it  was  in  the  one -dimensional  problem  previously  treated.  After  50  cn 
there  is  a  strong  departure  from  one-dimensional  behavior.  Figures  9(b)  and  9(c)  show 
the  development  of  transverse  reshaping  and  self  -  focus  sing .  The  pulse  characteristics 
as  a  function  of  propagation  distance  z  arc  summar  zed  in  Fig.  c'(d).  For  the  first  40 
centimeters,  the  area  and  intensity  along  the  axis  p  0  do  not  differ  significantly  lrom 
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their  behavior  in  the  one  dimensional  case.  The  half-width  in  the  p  dimension  is  de¬ 
creasing,  however,  indicating  transverse  reshaping  without  transverse  energy  flow. 

When  this  reshaping  process  has  reached  a  critical  point  there  is  a  sudden  increase  in 
the  area  and  energy  density  along  the  central  axis.  This  shows  that  there  must  be  sub¬ 
stantial  transverse  energy  flow.  The  energy  density  goes  through  a  peak  at  145  cm 
where  it  is  seven  times  its  input  value.  After  this  point,  it  starts  to  decay.  This  sub¬ 
sequent  energy  loss  is  not  accomplished  by  a  full  expansion  of  the  beam  size.  It  is 
noteworthy  that  substantial  ringing  has  occurred  after  the  pulse  peak  --  although  this 
effect  was  absent  in  the  one -dimensional  case.  The  energy  density  plotted  in  Fig.  )(d) 
is  integrated  up  to  20  nanoseconds  from  T  =  0,  and  thus  ib  an  underestimate  The 
existence  of  transverse  energy  flow  is  clearly  illustrated  in  Fig.  10  which  gives  the 
amplitude  and  phase  of  the  complex  field  amplitude  as  a  function  of  p  after  two  propaga¬ 
tion  distances,  z=  97  cm  and  z=  145cm.  For  each  propagation  distance  the  field  ampli¬ 
tude  and  phase  are  examined  for  two  values  of  T  ,  one  at  the  peak  ot  the  pulse  and  he 
other  shortly  after.  A  downward  curvature  to  the  phase  as  a  function  of  p  corresponds 
to  a  focussing  field.  At  the  97  cm  distance,  the  phase  variation  is  considerably  stronger 
after  the  peak  than  on  the  peak.  Had  the  remaining  propagation  been  in  free  soare,  the 
cut  through  the  peak  would  have  focused  after  675cm  whereas  the  cut  after  the  pe«k 
would  have  focused  after  only  107  cm.  The  actual  focus  occurred  after  40  cm. 

The  self-focussing  effects  just  described  disappear  when  the  transverse  relaxa¬ 
tion  time  is  shortened  to  equal  tne  pulse  duration,  keeping  other  parameters  constant. 
When  the  pulse  is  short  compared  to  the  transverse  relaxation  time  and  the  carrier 
frequency  is  sufficiently  offset  with  respect  to  tie  central  resonance  frequency  opposite 
signs  of  the  curvature  of  the  phase  as  a  function  of  p  develop  on  either  side  of  the 
resonance,  consistent  with  the  results  of  Grischhowsky  and  Armstrong. 

3 . 4  Summary 

We  have  studied  the  temporal  and  transverse  spatial  evolution  of  optical  pulses 
in  resonant  media.  N'unerical  integration  of  the  coupled  field-material  system  equations 
shows  the  phenomena  of  transverse  reshaping  and  self-fo'-ussing.  Self-focus  sing  is 
associated  with  transverse  energy  flow.  This  effect  is  specific  to  coherent  interactions 
between  the  pulse  and  resonant  medium,  and  disappears  when  the  pulse  duration  is 
comparable  to,  or  longer  than,  the  transverse  relaxation  time.  The  effect  is  distinct 
from  self-focussing  due  to  adiabatic  follvving  since  it  occurs  when  the  carrier  frequency 

is  at  line  center. 

The  implications  of  this  study  for  the  BALAD  receiver  are  discussed  in  Sec  lion 
4. 
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4.0  POSSIBLE  IMPROVEMENTS  IN  THE  DALAD  RECEIVER  SYSTEM 


4.  1  Comments  on  Alternative  BALAD  Geometries 

A  previous  report*  discussed  alternative  geometries  for  a  BALAD  receiver  sys¬ 
tem.  The  basic  features  of  the  BALAD  receiver  are  shown  in  Fig.  11(a),  and  a  rather 
impractical  multistage  absorber  alternative  is  shown  in  Fig.  11(b).  The  alternatives 
which  were  studied  in  particular,  were  the  optical  waveguide  configuration  and  the 
Fly's  Eye  lens  train  configuration  shown  in  Figs.  12  and  13.  In  both  of  these  cases,  the 
use  of  a  waveguiding  structure  in  the  saturable  absorber  allows  the  use  of  a  lmuh 
shorter  focal  length  lens  (i.e.,  a  smaller  f/number  lens)  at  the  input  to  the  saturable 
absorber.  Because  of  this,  the  signal  beam  cross-section  is  much  smaller  and  the 
signal  beam  intensity  is  correspondingly  larger  than  would  otherwise  be  the  case.  As 
a  result,  substantial  improvements  were  predicted  in  the  performance  of  a  BALAD 
system  employing  these  alternatives. 

The  present  contract  concentrated  on  the  demonstration  of  a  Xe-Xe  BALAD  sys¬ 
tem.  Thn  saturable  absorber  for  this  system  is  a  gaseous  discharge.  Because  of  the 
difficulty  in  obtaining  quiet  operating  conditions  in  the  xenon  absorber  discharge,  it  is 
apparent  that  the  optical  waveguide  geometry  would  not  be'appropriate  for  the  Xe-Xe 
BALAD  receiver  system.  It  would  be  difficult,  if  not  impossible,  to  run  a  discharge 
in  a  single  long  small  diameter  capillary  tube  of  the  type  .vhich  would  form  the  optical 
waveguides.  Based  on  our  present  experience,  such  a  discharge  (in  xenon)  would  sure¬ 
ly  be  noisy.  Finally,  the  array  of  waveguides,  which  would  require  running  many  such 
discharges  in  parallel,  would  present  an  impossible  task. 

A  passive  absorber,  such  as  the  SFb  employed  in  the  CO^-SF^  BALAD  system, 

is  clearly  compatible  with  the  "array  of  waveguides"  geometry.  However,  with  stub  a 

system,  one  would  be  tempted  to  provide  a  matching  array  of  waveguide  amplifiers  to 

precede  the  waveguide  absorber,  to  take  advantage  of  the  larger  gains  achievable  in 
1  1 0 

C02  wavei  uide  lasers.  If  so,  the  difficulty  of  running  parallel  discharges  would  have 
to  be  faced . 

The  Fly's  eye  lens  train  geometry  would  appear  to  be  t  ompatible  with  a  xenon 
discharge  absorber.  The  generally  observed  l/d  dependence  of  the  absorption  coefficient 
of  the  Xenon  absorber  must  be  recognized  in  the  design  of  such  a  system.  This  tends 
to  restrict  the  absorber  aperture  Finally,  the  task  of  producing  arrays  of  Fly's  eye 
lenses  using  infra-red  transmitting  materials  is  non-trivial. 
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Fig.  11(a)  -  BALAD  receiver.  Absorber  becomes 
transparent  at  the  focus  and  frequency 
of  the  signal  beam  . 


ABSORBER  CELLS  WITH  FIELD- 

/  .  \  \  I 


LFNS  WINDOWS 


LASER 

AMPLIFIER 


DETECTOR 


Fig.  11(b)  -  BALAD  receiver  with  fine-focus  multi¬ 
stage  absorber  required  for  gas  with  low 
attenuation  and  high  saturation  intensity. 
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4.2  The  Application  of  Non-linear  Sen-focussing  to  a  Geometrical  Loss  Filter 

As  discussed  in  Section  3,  the  analysis  of  the  propagation  of  a  coherent  optical 
pulse  in  a  resonant  absorber  revealed  a  strong  self -foe ussing  effect.  In  other  words, 
the  presence  of  the  optical  puls.'  locally  modified  the  focussing  properties  of  the  ab¬ 
sorber.  The  effect  was  strong  enough  to  focus  the  optical  pulse  into  small  spot  when 
otherwise  it  would  have  been  diverging. 

If  the  diameter  of  the  optical  beam  at  the  input  to  the  saturable  absorber  is  re¬ 
duced,  then  it  should  be  possible  to  find  a  beam  diameter  for  which  the  focussing  effect 
is  just  balanced  by  the  spreading  due  to  diffraction.  If  so,  then  the  nor -linear  self¬ 
focussing  effect  would  guide  the  optical  pulse  with  a  constant  diameter.  If  such  is  the 
case,  then  the  design  of  a  geometrical  filter  for  improving  the  performance  of  the 
BALAD  receiver  is  straightforward.  In  essence,  the  non-linear  self-focussing  effect 
would  provide  an  automatically  aligned  optical  waveguide  which  exists  only  for  the  sig¬ 
nal  beam  and  not  for  the  background  noise. 

Further  work  is  required  to  completely  analyze  the  behavior  of  such  non-linee r 
self-focussing  and  to  experimentally  test  the  propagation  of  optical  pulses  in  saturable 
absorbers.  The  purpose  of  the  following  discussion  is  to  anticipate  what  can  be  done  if 
that  work  were  successfully  completed  Thus,  the  discussion  may  be  useful  in  justify¬ 
ing  support  for  such  further  work. 

4.3  Design  of  a  Geometrical  Loss  Filter  Using  Non-linear  Self- foe  ussing 

The  essential  features  of  a  geometrical  loss  filler  using  non-linear  self-focus¬ 
sing  can  be  develo  >ed  from  the  BALAD  arrangement  shown  in  Fig,  14. 

Fi-st,  a  few  general  comments  about  the  overall  optical  system  are  in  order. 

For  simplicity  of  analysis,  the  figure  is  drawn  with  the  distances  L  =  L_  =  L  ,  and 

12  3 

the  lens  apertures  Dj  =  D^  =  D^  =  D^.  More  generally,  the  receiving  aperture  D 
would  be  larger  than  the  laser  amplifier  aperture  which  determines  D^  and  D^.  Also, 
the  length  L^  would  be  adjusted  in  order  for  the  beam  size  at  lens  No.  4  to  match  that 
beam  diameter  which  results  in  self-focused  waveguiding  in  the  saturable  absorber. 

It  is  clear  from  Fig.  14  that  the  desired  performance  of  the  indicated  optical 
system  is  as  follows:  a  signal  beam,  incoming  from  a  direction  0^  is  focused  to  a 
(diffraction  limited)  spot  at  lens  No.  2,  displaced  a  distance  X^  L_,  0^  from  the  optical 
axis.  Thus,  the  linear  field  of  view  of  the  receiver  is  given  by  =  D^/Lj.  If  lens 
No.  2  images  aperture  Dj  on  aperture  D^,  then  all  rays  which  pass  through  D(  (within 
the  field  of  view)  will  also  pass  through  D^.  The  spot  on  lens  No.  2  is  imaged  by  lens 
No.  3  at  lens  No.  4  where  it  is  displaced  by  X  =  (L./ L-,)X,  from  the  optical  axis. 
Finally,  lens  No.  4  renders  the  central  ray,  of  the  bundle  of  rays  within  the  signal 
beam,  parallel  to  the  optical  axis.  Thus  the  central  ray  would  still  be  displaced  by 
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at  the  output  of  the  saturable  absorber.  The  desired  focal  lengths  of  the  various 
lenses  are  determined  by:  f^  =  L^;  l/f^  =  l/L^  +  l/L^;  ^3  =  +  1/^3  anc*  ^4  =  ^  3 

Also,  the  (useful  lens  apertures  are  given  by:  Dj  =  L^/L^  •  D^,  =  L^/L^  •  D^. 

Note  that  the  optical  arrangement  in  Fig.  4  shows  the  signal  beam  diverging  in  the 
laser  amplifier  so  that  at  the  amplifier  output,  the  signal  beam  has  expanded  out  to  fill 
the  full  laser  amplifier  aperture.  If  the  signal  beam  diameter  at  the  amplifier  input  is 

QIo  /  2 

d,  then  if  the  amplifier  gain  e  >  (D^/d)  ,  this  arrangement  minimizes  saturation 
of  the  laser  amplifier  while  maintaining  the  maximum  field  of  view  allowed  by  the  laser 
amplifier  aspect  ratio. 


As  mentioned  previously,  the  purpose  of  lens  No.  4  is  to  render  the  central  ray 
of  the  signal  beam  parallel  to  the  optical  axis.  This  is  accomplished  by  the  require¬ 
ment  f^  =  L^.  Beyond  lens  No.  4,  in  the  absence  of  waveguiding,  the  signal  beam  would 
form  a  diverging  wavefront  as  shown  with  dotted  lines.  The  angle  of  this  cone  of  diverg¬ 
ing  rays  would  be  the  same  as  the  bundle  incident  on  lens  No.  4,  which  is  D^/L^  or 
(L^/Lj)  •  (Dj/Lj).  On  the  other  hand,  if  the  non-linear  s elf-focussing  of  the  saturable 
absorber  produces  waveguiding,  the  signal  beam  would  be  confined  to  the  narrow'  bundle 
of  rays  outlined  by  the  solid  lines  in  Fig.  14. 


Now  the  evolution  of  the  geometrical  loss  filter  begins.  H  there  were  an  aperture 
at  the  output  of  the  taturable  absorber,  and  if  this  aperture  were  properly  position^ 
and  its  diameter  matched  the  signal  beam  diameter,  then  in  the  second  case  above, 
(non-linear  self-focussing)  the  entire  signal  beam  would  pass  through  the  aperture  to 
the  detector.  Whereas,  in  the  absence  of  non-linear  self-focussing,  only  a  small 
fraction  of  the  signal  beam  would  pass  through  and  be  detected.  The  fraction  of  the 
total  power  received  in  the  absence  o'  self-focussing,  would  in  fact,  be  given  by 
l/(number  of  resolution  elemen  s  10  the  BALAD  receiver  field  of  view).  That  is,  stub 
a  (hypothetical)  filter  would,  in  tact,  introduce  an  attenuation  due  to  geometrical  loss 
for  any  signal  (or  noise)  which  is  \  recisely  that  factor  which  is  desired  in  order  to 
optimize  the  performance  of  the  B  3LAD  receiver.  (See  Ref.  19). 

Now,  of  course,  so  far  the  filter  is  hypothetical,  since  it  requires  the  automatic 
positioning  of  an  iris  or  diaphram  to  follow  the  signal  beam.  The  means  for  achieving 
tliis  result  is  to  introduce  a  matched  array  of  apertures  at  the  input  as  well  as  at  the 
output  and  at  selected  locations  down  the  length  of  the  saturable  absorber.  The  geomet¬ 
rical  considerations  for  such  a  "geometrical  filter"  are  illustrated  in  Fig.  15. 

Consider,  for  example,  a  close  packed  array  of  circular  apertures  with  centers 
separated  by  twice  the  diame'cr  of  the  apertures.  The  "unit  cell"  of  such  an  array  is 
an  equilateral  triangle  with  a  circle  centered  at  each  apex;  the  radius  of  the  ciri  It  is 
one  fourth  of  the  side  of  the  triangle.  The  average  transmission  of  such  an  array  (the 
fraction  of  the  area  which  is  open)  is  T  =  tt/B./T  =  0.226.  A  rectangular  array  with  the 
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c)  PLACEMENT  OF  APERTURE  ARRAYS 


Fig.  15  -  Aspects  of  geometrical  filters 


same  nearest  separation  would  have  T  =  ir/l6  =  .196.  Either  type  of  array  has  an  ac 
ceptably  large  transmission  for  the  use  contemplated  here. 


If  a  pair  of  such  arrays  were  closely  spaced  and  perfectly  aligned,  the  pair 
would  have  the  same  transmission  as  a  single  array.  Since  the  transmission  of  a  suc¬ 
cession  of  optical  eleme its  in  series  is  equal  to  the  product  of  the  individual  transmis¬ 
sions,  we  would  conclude  that  the  transmission  of  the  second  array  is  unity  (T  •  T  = 

Ti  implies  T2  =  when  il  immediately  follows  the  first  and  is  perfectly  aligned.  Con¬ 
versely,  if  the  second  array  were  laterally  displaced  by  an  aperture  diameter  in  any 
direction,  the  transmission  of  the  pair  of  arrays  would  be  zero.  From  this  we  conclude, 
that  the  transmission  of  the  second  filter  is  zero  if  it  immediately  follows  tne  first  and 
is  perfectly  misaligned  (Tj  •  T£  =  0  and  T }  -  0  -»  T  ^  =  0). 

Now,  if  the  pair  of  arrays  were  separated  longitudinally,  it  is  clear  that  the 
transmission  of  the  pair  of  arrays  would  depend  strongly  on  the  direction  of  the  incident 
light.  Collimated  light  incident  normally  (if  the  two  are  perfectly  aligned)  would  exhibit 
nearly  the  same  transmission,  reduced  only  by  the  diffraction  spreading  between  the 
two  arrays.  Clearly,  collimated  light  travelling  at  an  angle  0  from  the  normal  would 
be  (nearly)  perfectly  blocked  if  tan  0  =  D/L  where  D  is  the  aperture  diameter  and  L  is 
the  separation  between  the  two  arrays.  Again,  diffraction  spreading  of  the  light  would 
produce  som  transmission. 


But  then,  collimated  light  travelling  at  about  twice  that  angle  (tan  ©2  =  ZD/L) 
would  again  be  perfectly  transmitted  by  the  second  filter  --  if  the  azimuthal  direction 
of  the  light  path  produces  an  overlap  with  a  neighboring  aperture.  So,  one  is  forced  to 

introduce  at  least  a  third  array  of  apertures  to  discriminate  against  those  off  axis  light 
rays . 

Before  proceeding  to  investigate  the  minimum  number  and  optimum  placement 
of  such  arrays  of  apertures  or  "masks,  »  let  us  consider  the  behavior  of  t  ie  proposed 

improved  BALAD  system  making  use  of  non-linear  self-focus  sing  in  a  saturable  absorb- 
er 


First,  we  insist  that  an  input  s.gnal  beam,  whatever  its  incident  direction  (or 
lateral  displacement  at  the  input  to  the  saturable  absorber)  be  transmitted  by  the  filter. 
Ihus,  the  signal  beam  diameter  must  be  large  enough  so  that  at  least  one  aperture  is 
illuminated  by  (nearly)  the  full  intensity  of  the  signal  beam:.  The  worst  case,  of  course, 
is  to  have  the  signal  beam  precisely  centered  betv  fen  apertures  --  at  the  center  of  a 
"unit  cell"  of  the  aperture  array.  In  this  case,  3  (or  4)  apertures  would  be  illuminated, 
but  the  traction  of  the  incident  power  falling  on  one  aperture  would  be  l/( 3 .  3) 2  =  l/ll 
for  the  close  packed  array  or  l/(3.83)2=  l/l  4.  6  for  the  rectangular  array.  This  loss 
of  signal  power  can  be  compensated  by  increasing  the  laser  amplifier  gain.  An  addition¬ 
al  consideration  is  the  resulting  intensity  distribution.  Clearly,  the  beam  transmitted 
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by  an  aperture  would  not  be  a  gaussian.  It  would  be  closely  approximated  by  a  untiorm 
cylindrical  beam  which  is  sharply  truncated  at  the  radius  of  the  aperture.  The  behav¬ 
ior  of  such  a  beam  within  the  saturable  absorber  would  have  to  be  investigated.  More 
precisely,  it  will  be  necessary  to  investigate  the  performance  of  an  "aperture  wave¬ 
guide"  which  is  formed  by  the  sequence  of  aligned  apertures  down  the.  length  of  the 
absorber. 

It  is  clear  that  in  the  absence  of  self-focussing,  there  is  a  "bundle  of  rays"  in 
the  signal  beam  so  that  precise  alignment  of  the  axis  of  the  aperture  array  with  optical 
axis  would  be  unnecessary.  There  will  be  some  ray  within  the  bundle  which  would 
propagate  down  the  axis  of  the  array  sequence.  However,  the  non-linear  self-focus  sing 
will  presumably  confine  the  beam,  and  cause  it  to  propagate  as  a  constant  diameter 
beam  in  the  direction  determined  by  the  inciuent  direction.  Thus  it  will  be  necessary 
to  carefully  align  the  sequence  of  ape.tures  so  that  it  is  parallel  to  the  optical  axis, 
and  the  effect  of  this  misalignment  should  be  investigated  theoretically. 

One  method  of  distributing  the  arrays  of  apertures  down  the  length  is  to  ut  e  a 
geometric  progression.  If  there  are  a  total  of  N  2  apertures,  so  that  there  are  N  aper¬ 
tures  along  a  line  drawn  in  an  array,  then  we  can  number  these  sequentially  as 
I  =  1.2.3,  N.  The  problem  is  to  block  the  path  from  any  aperture  (I)  in  the  first 

array  to  any  other  aperture  (I+ J)  in  the  final  array.  The  first  additional  array  if 
o laced  midway  (at  L/2)  will  block  paths  from  I  to  I  +  ( 1  +  2M)  for  all  integer  value,  of 
M  but  the  paths  from  I  to  1+  2M  would  be  open.  The  second  additional  array  placed  at 
L/4  (or  3L/4)  will  block  half  of  these,  leaving  paths  from  I  to  I  +  4M  open.  The  third, 

U  (°r  7L/W)  Wl11  block  half  of  these  leaving  paths  from  I  to  I  +  8M  open.  The  kth 
additional  array,  placed  at  a  distance  of  L/2K  away  from  any  of  the  previous  arrays, 
will  leave  open  only  those  paths  from  I  to  I  +  (2KM).  Clearly,  when  2K  >  N,  all  paths 
except  I  =  I  will  be  blocked.  Thus  for  this  method  of  placement,  one  needs  (2+  log  N) 
masks,  including  the  mask  at  the  input  and  output  of  the  saturable  absorber.  2 

For  example,  5  masks  are  required  for  N  =  8;  6  masks  for  N=  16,  and  7  masks 
for  N=  32  (which  would  be  a  1024  element  array).  Thus,  the  number  of  masks  required 

is  not  excessive,  and  the  minimum  spacing  between  two  masks  is  L/N  which  is  also 
rea  onable. 

Finally,  it  can  be  mentioned  that  the  construction  of  such  arrays  of  apertures 
is  much  simpler  and  less  exacting  than  the  construction  of  an  array  of  "fly.s  eye  lenses" 
(or  an  array  of  waveguides).  Fewer  masks  should  be  required,  and  the  problem  of 

mounting  and  alignment  is  certainly  no  more  severe  than  the  corresponding  problems 
for  the  lenses. 
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4,4  The  Geometrical  Loss  Filter  -  A  Universal  Attenuator 

The  fact  that  the  geometric  loss  filter  is  universal  is  also  notable.  That  is,  it 
will  attenuate  any  and  all  incident  signals  and  noise,  except  that  signal  for  which  the 
saturable  absorber  exhibits  non-linear  self-focussing. 

The  absorption  of  the  saturable  absorber  will  contribute  additional  attenuation 
for  signals  which  coincide  in  frequency.  However,  one  major  point  of  the  geometrical 
loss  filter  with  non-linear  self-focussing  is  that  it  should  operate  and  produce  adequate 
attenuation  with  a  length  of  absorber  which  is  much  shorter  than  that  which  would  be 
required  if  one  was  relying  on  attenuation  alone.  That  is,  the  absorption  of  the  saturable 
absorber  will  probably  be  small  compared  to  the  geometric  loss. 

The  geometrical  loss  of  the  filter  can  easily  be  calculated  from  Fig.  14.  Let 
us  recall  that  each  aperture  of  the  array  represents,  optically,  a  re  olution  element  of 
the  BALAD  receiver,  and  that  the  total  area  of  the  array  represents  the  field  of  view 
of  the  receiver.  From  Fig.  14,  the  bundle  of  rays  incident  on  an  a  erture  at  the  input 
to  the  absorber  cell  will  have  the  (linear)  divergence  angle  D3/L3,  thus  at  the  end  of 
the  cell  it  will  spread  out  to  illuminate  an  area  ( D 3  •  L/L3)2  .  v  /  4  at  the  detector 
Since  the  total  area  is  <ir/4)D4  =  N2  •  (unit  cell  area),  and  only  rays  within  one  aperture 
at  the  output  will  reach  the  detector,  the  transmission  of  the  filter  will  be: 


Jarea  of  one  aperture)  /  T  ■  cell  area 
(illuminated  area)  TOT  area — 


TOT  area 
V ilium,  area 


T  ' 


D  L  2 
4  3 


V'  T^r, 


where  T  is  the  average  transmission  of  one  array,  i.e.,  area  of  aperture/area  of  a 
unit  cell.  In  the  simplified  geometry  of  Fig.  14,  the  latter  factor  is  unity,  and  the 
attenuation  is  T/N  .  Now  noise  uniformly  illuminates  the  first  mask  (at  the  input),  so 
that  noise  is  attenuated  by  an  additional  factor  of  T. 

In  comparison,  a  signal,  when  present,  will  be  focused  to  a  spot  whose  size  is 
by  design  somewhat  larger  than  the  area  of  one  unit  cell.  So  the  signal  passes  entirely 
through  one  aperture  (3  or  4  apertures)  and  suffers  an  attenuation  of  T  at  the  input. 
But  due  to  the  non-linear  self-focussing,  it  should  suffer  only  minimaMoss  thereafter. 
Thus,  the  signal-to-noise  ratio  would  be  improved  by  the  factor: 

’2  .  t 

( 4b) 


(S/N) 

'  out 


N‘ 


in 


TT 


i  .  s_+  t2  °4L3  2 

s  T<7t  +  T  \T7~L  / 

where  the  added  term  in  the  denominator  is  due  to  that  noise  which  is  spatially  (and 
spectrally)  coherent  with  the  signal  which  will  ''benefit"  from  the  no,  -linear  self¬ 
locussing  of  the  signal,  and  only  be  attenuated  by  the  signal  factor  T 

s  ' 


V 
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The  t  oherent 


noi,e  represents  1/N2  X  (B,)/(BT0T)  of  the  total  noise  power,  where  B  *  bandwidth 

Of  the  Signal  and  BTQT  =  total  noise  bandwidth  which  is  determined  by  the  laser  pre¬ 
amplifier. 

Compared  to  the  theore^cal  ideal  signal-to-noise  ratio  (S/N  =  P/h  v  B  ),  the 
signal-to-noise  ratio  would  be:  s 


(S/N), 


'ideal 


i  +  t^tot  n  2 

^  b^~  j  \  d'JTT  ) 


For  a  molecular  gas  laser  amplifier  such  as  CO.,  or  even  for  an  atomic  laser 
amplifier  such  as  xenon,  there  are  a  multitude  of  transitions  which  exhibit  gain  and  thus 
BTOt/Bs  Canbe  lar8e>  which  will  degrade  the  signal-to-noise  ratio.  However,  the 
output  of  the  geometrical  loss  filter,  which  is  collimated,  is  easily  adapted  to  the  in¬ 
sertion  of  a  Fabry-Perot  etalon  which  can  easily  filter  out  all  but  the  "working"  laser 

line.  In  such  a  case,  BTQT  =  Bandwidth  of  laser  lines  which  are  transmitted  by  the 
etalon. 

Note  that  in  this  case  we  have  not  been  concerned  with  the  overlap  of  the  absorber 
with  all  laser  lines  (although  the  absorption  would  always  help).  The  added  term  in  the 
denominator  of  Eq.  (49)  is  the  penalty  for  having  the  multitude  of  apertures. 

In  contrast,  in  the  "standard"  BALAD  arrangement,  where  the  absorber  must 
do  all  of  the  attenuation,  the  existence  of  a  non-ovrrlapped  laser  line  is  disastrous  -- 
the  corresponding  expression  for  this  case  is 

S/N  =  -  z  1  etot 

1  +  N  •  -IU1 
B 

s 

wfuch  is  worse  by  thy  factor  N2.  but  in  this  case,  of  course.  It  would  only  be  the 
sum  of  lines  which  do  not  overlap  the  absorber.  Moreover,  the  problem  of  filtering  is 
more  severe,  since  the  noise  input  (and  output  from  the  absorber  cell)  is  not  collimated. 

Thus,  in  addition  to  its  other  attractive  features,  .he  geometrical  loss  filter 
would  solve  the  problem  of  filtering  out  radiation  from  a  1  but  the  working  line  in  tin 
iO.  6m  COz  laser  anr  lifier  and  in  amplifiers  in  the  3-5u  spectral  range. 


5.0  SYSTEM  HAZARD  ANALYSIS 

The  diiectly  applicable  section  of  MIL-STD-882  is  paragraph  5.8.2.  1,  Prelim¬ 
inary  Hazard  Analysis.  According  to  this  paragraph,  a  preliminary  hazard  analysis 
shall  be  performed  as  the  initir’  analysis  task  during  the  acquisition  of  a  system.  This 
analysis  shall  be  a  comprehensive  qualitative  study.  Such  information  shall  be  used  in 
the  development  of  safety  criteria  to  be  imposed  in  performance  o,  design  specification* 
The  most  relevant  areas  of  those  listed  to  be  considered  are:  (f)  Effect  of  transient 
current,  electro-static  discharges,  electromagnetic  radiation,  and  ionizing  radiation 
to  or  by  the  system.  Design  of  critical  controls  to  prevent  inadvertent  activation  and 
employment  of  electrical  interlocks,  and  (g)  Use  of  pressure  vessels  and  associated 
plumbing,  fittings,  mountings,  and  hold-down  devices.  In  addition,  for  related  appli¬ 
cations,  items  (h),  (i),  (j),  (n),  (o),  (p),  (q),  (s),  and  (t)  should  be  considered. 

The  device  under  development  is  a  coherent  laser  receiver.  Under  normal 
proper  operation  of  the  device,  there  is  no  optical  output  whatsoever.  It  is  conceivable 
that  if  improperly  used,  the  laser  preamplifier  could  be  made  to  oscillate.  In  the 
present  system,  it  would  almost  require  a  deliberate  act  to  set  up  oscillation  conditions. 
But  as  a  general  rule,  the  alignment  of  optical  windows  or  other  partial  reflectors 
normal  to  the  optical  beam  axis  should  be  avoided.  A  surface  from  which  such  unwanted 
reflection  might  occur  is  the  optical  detector,  so  the  surface  of  the  detector  and  any 
associated  filters  should  be  deliberately  misaligned  to  avoid  the  possibility  ot  feedback. 
The  BALAD  receiver  is  particularly  protected  from  such  reflections,  since  the  gaseous 
bleachable  absorber  would  strongly  attenuate  the  reflected  beam  unless  the  frequency 
of  the  detected  signal  happens  to  coincide  with  the  saturable  absorber  line  center. 

Thus,  the  hazard  of  eye  damage  would  be  practically  ,  m-existent  in  the  operation 
of  the  BALAD  receiver.  In  contrast,  this  could  be  a  serious  problem  at  the  site  of  the 
associated  laser  transmitter. 

The  BALAD  receiver  system,  as  currently  realized,  does  involve  gaseous  dis¬ 
charges.  These  present  two  types  of  hazards:  high  voltage  power  supplies  and  vacuum 
systems  requiring  optical  windows. 

Utmost  caie  is  called  for  in  handling  the  high  voltages,  and  all  standard  pre¬ 
cautions  and  the  use  of  interlocks  are  strongly  recommended.  The  level  of  the  high 
voltages  required  can  often  be  reduced  by  paralleling  sections  of  discharges.  This  pro¬ 
cedure  is  strongly  recommended  since  a  reduction  in  the  length  of  discharge  sections 
is  also  desirable  from  the  point  of  view  of  achieving  quiet  d  sc  barge  c  onditions.  Of 

course,  all  such  parallel  discharge  sections  must  be  separately  ballasted  to  avoid 
instabilities. 
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The  vacuum  system  considerations  relate  more  to  reliability  rather  th  .n  safety. 
Because  of  the  small  size  and  the  geometry  of  the  discharge  system,  the  danger  from 
the  implosion  which  results  from  failure  of  the  vacuum  system  are  generally  minimal. 
The  use  of  glass  or  other  insulators  is  unavoidable  in  parts  of  the  discharge  system, 
but  the  system  ruggedness  and  reliability  can  be  improved  by  the  use  of  metal  and 
metal-ceramic  structures  where  possible.  The  optical  windows  are  obvious  weak  links 
but  the  use  of  O-ring  seals  on  the  one  hand  and  metal  to  ceramic  sealing  techniques  on 
the  other  have  resulted  in  rugged,  vacuum -tight  seals.  Alkali  halides,  which  are  often 
used  as  infra-red  optical  components  because  they  are  relatively  cheap,  are  to  be 
avaided  becauso  they  are  fragile. 

The  BAL^vD  receiver  is  a  precision  optical  device.  Thus,  the  same  techniques 
used  for  optical  telescopes  should  be  employed  for  mounting  and  supporting  this  system 
Greater  care  must  be  exercised  at  the  site  of  the  associated  laser  transmitter  --  which 
requires  far  more  isolation  from  vibration  and  acoustical  noise  sources.  The  safety 
considerations  are  similar  for  the  transmitter  and  BALAD  receiver,  except  that  the 
eye  hazard  must  be  a  major  concern  at  the  transmitter,  but  not  at  the  receiver. 


6.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  essential  features  of  an  xenon-xenon  BALAD  receiver  system  have  been 
demonstrated.  The  measured  properties  of  the  system  are  given  in  Table  I.  As  dis¬ 
cussed  in  Section  2,  this  demor  Uration  has  been  hampered  by  the  noise  in  the  xenon 
discharge. 

The  noise  performance  of  the  xenon  discharge  can  be  improved  by: 

1)  using  well-regulated,  low-ripple,  power  supplies,  and 

2)  using  shorter  discharges  which  are  well  ballasted. 

The  xenon  cells  can  be  sealed  off  from  the  vacuum  station,  and  the  xenon  pressure  can 
be  stabilized  by  controlling  the  temperature  of  a  liquid-solid  xenon  reservoir.^  The 
noise  performance  of  the  He-Xe  amplifier  at  3.5g.m  has  been  measured.  Under  proper 

.  i  1 

operating  conditions  the  effective  noise  per  mode  is  11  t  0  9x  10  '  which  is  close  to 

2 1 

ideal.  Thus  the  He-Xe  amplifier  is  a  good  quantum  amplifier  under  proper  operating 
conditions . 

Because  the  greater  difficulty  was  encountered  with  the  absorber  discharge,  a 
major  effort  should  be  expended  to  improve  that  discharge.  Clearly,  the  system  could 
be  improved  substantially  if  the  xenon  absorber  were  replaced  by  a  passive  absorber. 
Thus,  the  search  for  suitable  saturable  absorbers  at  3.  5u  and  at  other  wavelengths 
should  also  continue. 

The  theoretical  analysis  of  the  propagation  of  short  coherent  pulses  in  resonant 
absorbers  has  revealed  a  non-linear  self-focussing  effect.  This  effect  would  make 
possible  a  geometrical  loss  filter.  A  modified  BALAD  receiver  system  employing  sin  h 
a  filter  has  been  analyzed.  This  system  would  provide  substantial  direct  filtering  of 
the  radiation  from  other  transitions  in  the  CO^  or  other  laser  preamplifier,  and  would 
facilitate  the  use  of  standard  filters  for  eliminating  the  remainder  of  such  broadband 
background  radiation.  Thus,  the  further  study  of  non-linear  self-focussing  is  recom¬ 
mended,  to  be  followed  by  the  development  of  the  geometrical  loss  filter  if  the  result- 
of  the  study  are  favorable. 
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TABLE  I 


Xe-Xe  BALAD  RECEIVER  CHARACTERISTICS 


Input  signal  from  a  He-Xe  laser  at  3.  5pm 

Output  signal  to  the  detector 

Amplifier  length 

Amplifier  diameter 

Absorber  length 

Absorber  diameter 

Gas  pressure  in  He-Xe  amplifier 

Gas  pressure  in  Xe  absorber 

Small  signal  gain  coefficient  of  amplifier 

Small  signal  absorption  coefficient  of 
absorber  a 

o 

Saturation  intensity  of  amplifier 
Saturation  intensity  of  absorber 
Amplification  bandwidth  at  3.  5pm 
Absorption  bandwidth  at  3.  5pm  (A\^) 

Saturated  transmission  hole  width  in 
absorber  B 

Possible  number  of  spectral  modes  (AvD/B) 

Figure-of-Merit  of  absorber  [F=(q  B/l  )] 

O  '  S/J 

Possible  number  of  spatial  modes  (resolution 
elements)  of  amplifier 

Possible  number  of  spatial  modes  (resolution 
elements)  of  absorber 

Minimum  spot  size  in  absorber 


0.  28.  W 
6.7  pW 
68  cm 
1  cm 
20  cm 
0.  5  cm 

2  Torr  He,  0. 1  Torr  Xe 
0.  1  Torr  X  e 
-JldB/m 

-60  dB/ m 

-2.5  mW/cm2 
>0.  51  mW/cm2 
*  1 00  MHz 
- 1  00  MHz 

*■  1  6  MHz 
6 

2.  9  x  1  0*  m/joule 

■4H 

7 

0.17  mm 
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Appendix  A 

Propagation  of  a  C  W  Gaussian  Beam  in  a  Saturable  Absorber 


The  propagation  of  a  continuous  wave  ("c\v")  confocal  beam  which 
has  a  Gaussian  radial  profile  has  been  studis  d  to  investigate  the  possibility  of 
self-focusing  of  such  beams.  This  analysis  is  based  on  the  formalism  which 
was  developed  to  describe  the  propagation  of  a  non- saturating  beam  in  a  medium 
with  a  specified  radial  dependence  of  gain  coefficient  a  and/or  index  of  re¬ 
fraction  n.  In  order  to  present  the  formalism,  a  brief  review  of  the  latter 
problem,  will  first  be  given.  A  parabolic  radial  variation  is  assumed  for 
the  gain  coefficient  and  index  of  refraction 


a  =  a 

o 


2  a2r 


(A.  1) 


n  =  n 

o 


1  2 
2  n2r 


(A.  2) 


where  the  index  of  refraction  is  proportional  to  the  (real)  propagation  constant  c 


^  '  So  "  iS2r2  =  (2rT/X)n  '  (A.  3) 

For  an  abs^j-her,  a  is  negative.  L.  Casperson,  et  al/2^  extended  H.  Kogel- 
nik's  result  to  include  the  frequency  dependence  of  both  gain  and  index  of  re¬ 
fraction  (that  is,  the  "anomalous  dispersion")  to  apply  to  a  high  gain  3.  5  um 
xenon  laser  operated  near  threshold.  A.  Javan,  et  al.(j)  studied  the  problem 
of  self- focusing  due  to  intensity  dependent  anomalous  dispersion.  In  their 
treatment  of  the  problem,  the  inhomogeneous  broadening  was  approximated  by 
a  Lorentzian.  In  the  following,  the  propagation  of  a  cw  Gaussian  beam  is 
studied,  taking  into  account  both  a  specified  radial  dependence  of  gain  or  ab¬ 
sorption  and  index  of  refraction  such  as  that  caused  by  a  discharge  and  also 
the  radial  dependence  caused  by  saturation  by  the  Gaussian  beam  profile. 

It  is  convenient  to  combine  (A.  1)  and  (A.  3)  by  using  the  complex 
propagation  constant  k  ?  4-  ia  ,  i.  e.  , 


where 


A-l 


k  =  S  -t-ift  and  k-  =  B.+ia. 
o  o  o  2  2  2 

Assuming  the  Gaussian  beam  is  of  the  form 

l  k  7 

E  =  cxp{-i[koz  +  t(z)]}  exp[-i  ^  ;^yr  ] 


(A.  5) 


where 


kQ  is  assumed  to  be  a  constant  independent  of  z, 
tjr(z)  is  *  complex  phase  parameter, 


and 


n  2* 

=  yx  +y  is  the  distance  from  the  optic  axis, 


q(z)  is  the  complex  beam  parameter  given  by 


1  .  \ 

-  l 


where 


R<*>  ™2<z, 


(A.  ..) 


and 


R(z)  is  the  radius  of  curvature  of  the  beam, 
UJ ( z )  is  the  spot  size  of  the  beam, 

X  is  the  wavelength  of  the  beam. 


If  ti.e  Gaussian  beam  given  by  (A.  5)  is  substituted  into  the  scalar  wave  equation 

(A.  7) 


V2E  +  k2E  =  0 


wh  ere 


2  2  2 
k^  *  -  k  k.r 

o  o  2 


one  gets  the  Ricatti  equation  which  governs  the  propagation  of  the  complex  beam 
parameter,  i.  e.  , 


(~  )  I  -  0 

'  q  dz  q  k 

o 


(A.  8) 


and  the  relation  between  the  complex  phase  parameter  and  the  complex  beam 
parameter  is  given  by 


d  ,  ,  .  _  .1 

dz  ^  Z  '  "  1  q(z ) 


( A.  (>) 


For  a  given  input  beam  parameter  qlf  the  output  beam  parameter  can  be 
obtained  from  the  solution  of  (A.  8),  i.  e. , 


~~  cos^ik-  'k  ) z  -  Jk  /k  sin  7(k_/k  )z 
1  qi  2  °  2  o  2  o' 


lZ  1 


N/K7kT  sin  ./(k,  /k  Iz  +  cos  ,/(k_  /k  I 

CO 


(A.  10) 


q,  v  o  2 


Consider  the  following  cases: 

(a)  Jk^Tk~  =  ±  (6  +  ie ),  6  >  0,  e  >  0,  i.  e.  ,  6  and  e  have  the  same  signs 

where  6  and  e  are  the  real  and  imaginary  part  of  ^^/k  respectively.  Writing 
cosine  and  sine  functions  in  terms  of  exponential  function  in  (A.  10)  and  taking  the 
limit  z  -  ®,  one  gets  the  stationary  solution 


( ~)  =  -i(6  +  ie)  =  e  -  i6 


l2  z~>  ® 


(A.  1  1 ) 


From  the  definition  of  q  given  by  (A.  6),  the  stationary  radius  of  curvature  R£ 
and  the  spot  size  a1  are  given  by 


a  =  /  X  /tt  6 


(A.  12) 


(A.  1  > ) 


(h)  Jk^/kQ  =  ±  (6  +  ie  ),  6  >  0,  e  <  0,  i.  e.  ,  6  and  c  ha 
Similar  to  case  (a),  one  gets  the  stationary  value 

(~“)  =  i(6  +  ie )  =  -e  +  i6 

"2  z  -  ® 


ve  opposite  signs. 


(A.  14) 


(A.  15) 


U)  =  ,/X/-n6 

00 


(A.  I,  ] 


One  can  see  that  the  spot  size  uco  as  given  by  (A.  lo)  will  be  imaginary. 

(c)  ^2^0  ~  e’  '  e^^er  0)  a2  ao  e’  ’  no  Pa^n  or  absorption 

or  (2)  C^/CLq  ~  ®2^®o’  **  c*  ’  an  identical  radial  dependence  of  gam  or  absorption 
and  index  of  refraction.  (Note  that  2nn/X  i  0  always.  ) 


From  solution  (A.  10),  one  obtain,  that  for  any  input  beam  parameter 

the  output  beam  parameter  ,2  oscillate,  periodically  without  approaching  a 
stationary  value. 

(d)  7(kJ7k^T  -  0,  i.  e.  ,  (*2  =  0,  82  =  0,  no  radial  dependence  of  either 
gain  or  absorption  coefficient  and  index  of  refraction. 

From  solution  (A.  10).  one  gets  the  free  space  result,  i.  e.  . 


q2  =  ql+  2 


Consider  the  expression 


(A.  17) 


A  - 

V  k 


/B2  +  ia2 
9o+lao 


1  +  i(a  /8  )  / 

o  o 


where 


*2,2  l'4 

-  (-  l  2  N  i* 

e2  +a2  e 

o  o 


i  1/,  -1  ®2  -1  ®2  i  CC  v 

6"  +tan  12  -  tan"1  ^  ) 
o  °2  Bo  ' 

Hence  the  real  and  imaginary  part  of  can  be  written  as 


g  2  ,  2  1/4 

o  o 


(A.  18) 


(A.  10) 


e2  2  [/4 

-  (  2+aZ  \ 

~  V7?"j 

o  o 


It  can  be  shown  that  the  condition  for  6  >0  and  c  >  0  is  given  by 


a2>fl2^ 


I  A,  201 


(A.  21) 


Rad,al  dependence  of  gain  or  absorption  coefficient  and  index  of  refraction 
can  be  caused  by  either  the  discharge  in  the  plasma  tube  or  the  saturation  bv 
a  radial  dependence  Caussian  beam  or  the  combination  of  the  above  two  effects. 


A-4 


The  radial  dependence  of  gain  coefficient  due  to  discharge  can  be  assumed  to 
have  the  form^ 


a(x)  =  a 

o 


(A.  22) 


where 


a  (x)  =  aDe 
o  P 


-x 


V  -  V 

x  =  2jlnZ 


Av 


D 


(A.  23) 


Av 


D: 


*  o' 

In  the  case  of  absorber, 
have 


the  peak  gain  value  at  line  center 

the  full  doppler  width  at  half  the  central  val.ie 

the  radius  of  the  plasma  tube. 

ap  is  negative.  Comparing  (A.  22)  with  (A.  1),  we 


fto  =  ao(x) 

2ao(x) 

a2  *  2 
r 

o 


(A.  24) 
(A.  25) 


The  corresponding  index  of  refraction  can  be  calculated  by  Kramers  Kronig 
relations 


(n  -  1 )  =  -^  P.  V. 

A.  TT 


tt(v') 


V  -  V 


dv' 


2ap  /  2 


where  F(x)  is  Dawson’s  integral^  ^  '  given  by 


F(x)  =  e 


•x 


2  x  2 
et  dt 


The  index  of  refraction  can  be  written  as 


(A.  26) 


(A.  27) 


(A.  28) 


iTH 
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n(x)  =  1  + 


(A.  29) 


Comparing  (A.  29)  with  {A. 2  ),  we  have 


n 

o 


1  f 


Xap 

~m  F(x) 


n2  = 


2>ap 

TTiT 


F(x) 


(A.  30) 


(A.  31) 


The  radial  dependence  of  gain  or  absorption  and  index  of  refraction  caused  by 
saturation  of  a  Gaussian  beam  can  be  calculated  as  follows.  The  gain  coeffi¬ 
cient  of  the  doppler  broadened  medium  in  the  presence  of  a  saturation  field  I '  is 
given  by 


a(x)  = 


ao(x) 


yi  +  (i7is) 


(A.  32) 


where  broad  doppler  width  has  been  assumed.  I  is  the  saturation  intensity  of 

s  7 

the  medium.  Let  the  Gaussian  b»am  be 


I'  =  Iexp[-(r2/u:2)] 


(A.  33) 


where  U)  is  the  spot  size  which  reduces  the  intensity  by  1/e  of  its  central 
value.  Substituting  (A.  33)  into  (A.  32)  and  making  an  expansion  in  the  region 
(r2/UJ2)«  1,  one  gets 


a(x) 


goW  r 
v’l  +  (I/Is)  L 


1 


1  1/1 
—  s 

2  1  +  (I/Is) 


Comparing  (A.  34)  with  (A.  1  ),  the  following  results  are  obtained 

a  ao(x) 

a°  "  71+ (I/I 

s 

ao(x)  i/lg  ! 

“2  =  'yTT7T7rT  TtW7  ' 

s 


(A.  34) 


(A.  33) 


/A.  3(, ) 
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The  index  of  refraction  for  inhomogeneous  broadening  laser  medium  in  the  limit 


2^(-£r)'/TT^7y  «  1  , 


as 


where  Av  is  the  full  homogeneous  line  width  at  half  the  central  value  can  be  written 

(5) 

n(x)  =  1  .  2  ,/u  xexp(-x2)]  .  (A.  37) 

.At  d 


Substituting  the  Gaussian  beam  (A.  33)  into  37),  we  get  the  following  ex- 

2  2 

pression  in  the  limit  (r  /w  )«  1, 

•'(x)  =  1  2-/<^(^'Vlt(1/1B>  (UI  i  i  (I/Ia)  ^-)*exp(-x2)]. 

(A.  38) 


Comparing  (A.  38)  with  (A.  2  ),  one  gets 

Xa 

n0  ~  1  + -  2  ^  (m^) '/TTW  x^P<-x2)]  <A-«) 


Xa 


'D 
I/I 


n2  = 


_P_  2  (-£*-) 

'  Avc  1  +  (I/I  ) 

s 


s _ x  exp(-x  ) 

2 

0) 


(A.  40) 


If  radial  dependence  of  gain  or  absorption  and  index  of  refraction  clue  to  both 
discharge  and  the  saturation  effect  are  taken  into  considerations,  then  the  gain 
coefficient  and  the  index  of  refraction  are  given  by 


a  = 

o 


ao(x) 


a2  =  ao(x)Pf - 


n/1  +  (I/Is) 
I/I 


r  2  [1  +  (I/I  )]3/2  u,2 


) 


(A.  4  1 


(A.  42) 


Xa, 


no  =  1  2'/Ji:2'l-Srr)'/rnI7V  xeXp(-x2)]  (A.  4 


4  3) 


D 


Xa, 


I/I 


n2 


C_2FlxL  .  - JL_  afc  )  1 

lvffir2  VlvD  ,/l  i  U/l  )  U)2 


(A.  i  n 


In  the  following,  the  focusing  and  defocusing  effects  of  saturable  absorber  will 
be  studied.  The  results  will  be  compared  with  those  obtained  from  an  amplifier. 
Generally,  8o  is  much  larger  than  aQ.  i.  c.  ,  8o  7t>  a  ,  one  can  write 


,/k^rr  *  V82  +  ia2/8o  =  jqj r  yrrr-i(62/a2)] 


(A.  45) 


The  ratio  82/a2  can  be  shown  to  be  much  smaller  than  one  for  small  frequency  devia¬ 
tion  from  the  line  center  for  the  following  two  cases: 

Case  1:  and  6^  caused  by  the  discharge. 

Taking  the  ratio  of  (A.  31)  and  (A.  25),  we  have 


9  2  (2n/>ln  2  2  FM 


2.  2  K/n  exp(-x  ) 

(A.  46)  is  smaller  than  one  if  x  is  much  smaller  than  about  1/2. 

Case  2:  and  P2  caused  by  saturation  of  a  Gaussian  beam. 

From  (A.  36)  and  (A.  40),  we  have 


(A.  46) 


(A.  47) 


(A.  47)  is  much  smaller  than  one  if 


[1  +  (I/I  )]  4  v rin2  (Av/Av_.) 

s  L) 

^or  small  frequency  deviation  from  the  line  center,  expression  (A.  45)  can  be 
expanded  in  series. 

Consider  the  case  where  a  weak  Gaussian  beam  is  propagating  through  an 
absorber  which  has  the  absorption  coefficient  and  index  of  refraction  given  by 
(A.  22)  and  (A.  24)  respectively.  In  this  case,  ap  given  by  (A.  23)  is  negative. 

From  (A.  24)  and  (A.  25),  the  absorption  coefficient  of  the  absorber  has  the 
following  sign: 

aQ  <  0  and  <  0  . 
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(i)  For  small  frequency  shift  to  the  high  frequency  side  from  the  atomic 
line  center,  i.  e. ,  x  >  0,  the  radial  dependent  par  of  refractive  index  n  has 
a  positive  sign  as  can  be  seen  from  (A.  31).  This  implies  that  >  0*  Writing 
a_  =  -  |  aJ  ,  we  have  from  (A.  45)  the  following  result 

fa  fa 

./iyiT  .  /l^TTir  ^  (i  +  i  i  -^y) 

Expression  (A.  48)  is  valid  for  02/|a2|<<*l  • 

The  real  and  imaginary  part  of  <fk are  given  by 

6=/T5772T  (ui^) 

' =  -  Vlo2l/zeo  (‘--Jijg-)  • 

6  and  e  have  opposite  signs.  If  (A.  49)  is  substituted  into  (A.  16),  one  sees  that 
ai^  will  be  imaginary.  Focusing  cannot  occur. 

(ii)  For  small  frequency  shift  to  the  low  frequency  side  from  the  line  center, 

9  will  be  negative,  i.  e.  ,  9  <  0.  Writing  9  =-|8_|,  we  have  from  (A.  45 )  in  the 

,.2  .  |P2|  2  Z  Z 

limit  I  c  i.-  «  l  f 

^  yk2/ko  =V|a2|/28o  ‘(1  +  i  T^j)]  •  (A-5" 

Therefore  6  and  e  are  given  by 

|  q  | 

6  =  •/  I  a2  |  /28o  (l  -  “  y^)  (A.  52) 

€  -  -  V|a2|/2«o  +  \  T^]-)  •  (A.  5i) 

Again,  f  and  e  have  opposite  signsso  that  focusing  cannot  be  obtained. 

Consider  the  case  where  the  absorber  is  saturated  by  a  strong 
Gaussian  beam.  The  radial  dependence  of  absorption  coefficient  and  index  of 


-ilk)]- 


(A.  48) 


(A.  49) 


(A.  50) 
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refraction  are  assumed  to  be  mainly  produced  by  the  strong  field  so  that  (A.  34) 
and  (A.  38)  are  applicable.  The  absorption  coefficient  of  the  saturable  absorber 
have  the  following  signs:  aQ  <0,  0L,  >  0.  aQ  and  a.,  are  given  by  (A.  35)  and  (A.  36), 
respectively,  where  G,q(x)  is  negative. 

(iii)  On  the  high  frequency  side  from  the  line  center  >  0),  we  can  obtain 
from  (A.  45)  in  the  limit  <<:"*  ’ 


6  =  ,'a2/29o  [l  +  (*2/a2)J 

e  =  ./a2/2Sori-(e2/a2)J 


6  and  €  have  the  same  sign.  There  is  a  stationary  value  of  Focusing  will 

occur.  From  (A.  12),  (A.  13),  (A.  54),  and  (A.  55).  one  can  obtain  the  expressions 
of  stationary  spot  size  and  radius  of  curvature 


0) 


r  2 

(rr/Xr 


S2.2  -W4 

it  vl  +h:  '  J 

o  2 


l2  f , 


(A.  5o! 


(A.  57) 


If  the  expressions  of  aQ,  nQ  and  ?z/az  as  give-  by  (A.  3b),  (A.  30),  and  (A.  47 
arc  substituted  into  (A.  56)  and  (A.  57),  one  obtains  after  some  calculations 
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R  (x)  = 


4  j 1  +  _Z  -  2  JUH  Jl  +  \l/ls)  X  exp(-x2 )]  J 


Otp  exP("x  /  I/I  s 

IrriurT  ttW 


(A.  59) 


(A.  58)  and  (A.  59)  are  valid  for 

0  <  x  « 


(i  t  j-)  4 \/Tn2  (- 


(iv)  On  the  low  frequency  ^ide  from  the  line  center,  3^  <  0,  we  can  obtain 
from  (A.  45)  in  the  limit  ■  <<  1  , 

l°2l 


5  =  vW2^: 


“2  ' 

(A.  60) 

J**1) 

a,  / 

(A.  61) 

6  and  c  have  the  same  sign  so  that  focusing  occurs.  As  the  frequency  shifts 
more  to  the  lower  frequency  side  from  the  line  center,  |  |  will  increase  and 

6  will  decrease,  as  can  be  seen  from  (A.  bO).  Since  the  spot  size  varies  in¬ 
versely  as  the  square  root  of  6,  it  will  increase  as  5  decreases.  Similarly, 
the  stationary  spot  size  is  found  to  be 


uu  (x)  - 

CD  ' 


4  1  +  _^P[£^I  .  2,Tn2  (-7~^)  •/1  +  (1/lg)  N  exp(-x2lj  | 

^/apexp(-x2)  1/  lg  ^  r  /  T  N  /  £V  V.T 

tm^tT  ^  i. 


.  (A.  62) 


R  (x)  has  identical  express. on  as  (A.  59).  (A.  o2)  is  valid  for 
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(A.  58)  and  (A.  62)  can  be  combined  to  one  expression 


A  plot  of  (A.  63)  and  (A.  59)  are  shown  in  Figs.  A.  1  and  A.  2  for  the  case  of  xenon 
absorber  at  different  intensities.  The  following  value  of  parameters  are  used: 

X  =  3.  5am,  ap  =  13.  85/m,  Av/AvD  «  0.  1. 

As  a  comparison  to  the  amplifier,  one  can  easily  show  that  for  the  case 
of  unsaturated  amplifier  which  has  the  radial  dependence  of  gain  coefficient 
and  index  of  refraction  due  to  the  discharge  alone  has  the  same  expressions 
as  (A.  54)  and  (A.  55)  for  >0,  &2  >  0,  82  >  0  (i.  e.  ,  high  frequency  side  and 
the  same  expression  as  (A.  60)  and  (A.  61)  for  0q  >  0,  &2>  0,  ®  e’  • 

low  frequency  side).  Therefore,  focusing  can  be  obtained  on  the  high  frequency 
side  from  line  center  and  defocusing  will  be  obtained  on  the  low  frequency  side 
from  the  line  center.  The  above  results  have  been  pointed  out  by  L.  Caspersion, 
e,al.<2  » 

In  the  case  of  a  saturated  amplifier  whose  radial  dependence  of  gain  coef¬ 
ficient  and  index  of  refraction  are  dominated  by  the  saturating  Gaussian  field, 
identical  results  as  given  by  (A  49)  and  (A.  50)  will  be  obtained  for  3^  >  0, 

<  0,  3  2  >  °  frequency  side),  and  identical  results  as  given  by  (A.  ~>Z) 

and  (A.  53)  will  be  obtained  for  aQ  >0,  a2  <  0,  °2  <  0  (low  frequency  side). 
Therefore  focusing  will  not  occur  on  either  side  of  the  line  center. 

One  interesting  case  is  when  there  is  only  radial  dependenc'’  of  refractive 
index  and  neither  radial  dependence  of  gain  nor  absorption,  i.  e.  ,  5^  t  0,  3, 


A  -  l  2 


-0.8 

Fig.  A.  1 


i  I 

\ 


-0.6  -0.4  -0.2  0  0.2  0.4  0.6  0  8 

-  (Solid  Line)  Stationary  beam  radius  vs.  frequency  deviation  from  line 

center  (x  =  2./I3£T  for  l/lg  (=  0.  1,  l,  3),  P2  /  0  (i.e.  ,  radial 

dependence  of  index  of  refraction).  (Broken  Line)  i/l  =1,3-0 

'  s  2 

(i.e.  ,  no  radial  dependence  of  index  of  refraction). 
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mTN^lx) 


(Solid  Line)  Radius  of  curvature  of  the  Gaussian  wave  vs.  frequency 
deviation  from  line  center  for  i/l  (=0.  1,  1,  3),  P-,  4  0  (i.e.  ,  radial 
dependence  of  index  of  refraction).  (Broken  Line)  i/l s  =  1,  P-,  »=  0 
(i.e.,  no  radial  dependence  of  index  of  refraction). 


the  medium  behaves  like  a  converging  lens  only  on  one  frequency  side  from 
line  center.  In  the  case  of  absorber  (aQ  <  0),  there  will  be  focusing  on  high 
frequency  side  from  the  line  center  and  defocusing  on  low  frequency  side  from 
the  line  center.  The  stationary  spot  size  and  radius  of  curvature  for  a  <  0, 

>  0  (high  frequency  side)  can  be  found  to  be 


=  f(TT/X)2  (°2/P0)]‘1/4 

i  <lao|/Po)2  (S2/9o)] 


-1/2 


(A.  64) 


(A.  65) 


In  the  case  of  the  amplifier  (aQ  >  0),  there  will  be  focusing  on  low  frequency- 
side  from  the  line  center  and  there  will  be  defocusing  on  high  frequency  side 
from  the  line  center.  The  stationary  spot  size  and  radius  of  curvature  for 
ao  >  °'  <  °  ^low  frecluency  side)  can  be  found  to  be 

*.'[<"A)2](ao/0o)2(|‘(2|/9o)]  1/4  (A.  66) 


K  =  (te2|/20  )-|/2  . 


(A.  67) 


(3  ) 


The  above  properties  have  been  pointed  out  by  A.  Javan,  et  al. 1 

The  other  interesting  case  is  when  there  is  only  radial  dependence  of 
either  gain  or  absorption  and  no  radial  dependence  of  refractive  index,  i.  e.  , 
a2  *  °'  32  =  °'  the  medium  behaves  as  a  converging  lens  only  when  the  ampli¬ 
fier  is  not  saturated  and  also  when  the  absorber  is  saturated.  The  stationary 
spot  size  and  radius  of  cc 
a,  >  0  can  be  found  to  be 


spot  size  and  radius  of  curvature  of  an  unsaturated  amplifier,  i.  e.  ,  a  >  0, 


x»  f(^A)2(a2/2Po)r1/4 


R«  =  <a2/23o> 


-1/2 


(A.  -  8) 

(A.  or ) 


Results  (A.  08 )  and  (A.  6d)  have  been  calculated  by  H.  Kogelnik.*  1  ^  For  a 
saturated  absorber,  i.  e.  ,  aQ  <  0,  a2  >  0,  the  stationary  spot  size  and  radius 
of  curvature  are  identical  to  (A.  68)  and  (A.  69).  A  plot  of  (A.  68)  and  (A.  69) 
are  shown  in  Figs.  A.  1  and  A.  2  by  a  dash  line  as  a  comparison  for  the  cases  where 


=  2  /°. 
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Appendix  B 

Analysis  and  Description  of  Numerical  Procedure 

We  will  now  discuss  the  numerical  integration  of  the  normalized 
working  equations: 

a  i  =  ic.fa  +  ;bje-"c  +g^  (b.  i) 

a  9  =  -  {—  +  Auu)0  +  e  VV  (B.  2a) 

T  2 

a  VV  =  -(W-  W  )/T.  -  Re{^f  ’  |  (B.  2b) 

T  O  i 

where  c,,  5,  g,  T.,  T_,  Au),  VV  are  constants  and  9  and  £  are  the  normalized 
slowly  varying  complex  envelopes  of  the  macroscopic  polarization  and  electric 
field  respectively.  W  is  the  macroscopic  population  inversion  of  the  two-  state  atomic 
system.  The  numerical  scheme  we  have  chosen  is  basically  described  in  the  paper  by 
Icsevtii  and  Lamb.*  The  normalized  transformed  independent  variables  -  I  arc 

defined  in  the  text.  The  input  conditions  are  £  (z-0,  o  1  and  9,  W(’',  0,  »  -  0).  That  is. 
the  £ -field  is  prescribed  just  inside  the  beginning  of  the  medium  (the  wave 
propagates  i.i  the  positive  z-direction)  and  W,9  assume  their  initial  status 
before  the  arrival  of  the  perturbing  £-field.  Computation  always  begins  at 
-  -  0  which  physically  would  correspond  to  an  observer  at  some  position  z 
down  range  turning  on  his  detector  only  after  the  a  time  z/c  has  past  since  the 
pulse  entered  the  medium  at  z-0.  The  3-dimensional  grid  space  o  ,  r  ,  . 

t  where  the  field  is  to  be  evaluated  is  filled  in  a  manner  consistent  with  tne 

n 

proper  use  of  a  predictor-corrector  numerical  integration  scheme  for  both 
the  field  equation  and  medium  equations.  Once  the  £ -field  is  known  at  some 
z-plane,  we  use  Eq.  (R.  1)  to  guess  (predict)  the  £-field  in  the  entire  z  -  JLz 
plane.  The  value  of  Az  is  necessarily  small  since  we  are  assuming  an  i  n 
variations  are  small  between  z  and  z  t  fz.  Now,  given  the  values  of  9  and 
VV  at  ”  0  for  all  P,  we  march  9  and  w  in  r  using  Eqs.  (B.  2a)  and  (B.  2b)  as 

if  they  were  initial  value  ordinary  differential  equations.  The  values  of  the 
mesh  separations  are  determined  by  the  combined  considerations  of  local 
stability  tests,  experimentation,  and  cost  effective  use  of  computation  lime. 

The  dependent  variable  guessed  (predicted)  at  the  next  forwai  !  mesh  point  is 
evaluated  usinfr  only  quantities  that  have  been  corrected.  For  example,  to 
predict  9  at  some  time  point  t  il’  only  previously  corrected  values  of 
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9,  l,  and  W  at  mesh  points  t  and  T  -  At  are  used.  This  is  absolutely  necessary 
if  the  predictor-corrector  (P-C)  is  to  be  used  at  maximum  accuracy.  We  have 
analyzed  the  scheme  and  found  it  to  be  conditionally  stable  and  consistent  to 
within  2nd  order  in  At,  Ao,  A^.  A  more  detailed  derivation  of  these  results 
is  given  in  N.  Wright's  Fh.  D.  thesis  2  .  It  is  known  for  maximum  ef¬ 
ficiency  of  the  P-C  scheme  one  should  use  only  one  iteration.  This  we  have 
done.  We  found  only  one  special  case  where  the  9  and  W  equations  become 
unstable  (namely,  the  rate  equation  limit  of  very  short  T^Tp/  5).  The  st  ability  prob¬ 
lem  always  arose  in  the  C-field  equation  since  the  ratio  of  c  )/(Z?)  does 
have  values  where  the  scheme  is  conditionally  unstable.  For  maximum  ac¬ 
curacy  the  parameter  A*/Ao2  was  chosen  to  be  as  close  to  the  unstable  region 
as  possible.  Other  researchers  have  also  found  this  condition  to  be  true. 

We  have  used  the  well-known  trick  of  transforming  to  a  stretched  coordinate 
framed  *  Stability  considerations  are  uneffected  by  such  a  transformation  since 
the  stability  is  evaluated  locally  where  the  equations  are  linear.  The  trans¬ 
formation  we  used  clustered  sampling  data  where  the  derivatives  with  respect 
to  t  were  largest.  The  explicit  transformation  was:  T  =  1/2(1  +  tanh  { 0 
where  a  is  the  stretching,  factor  making  mesh  points  more  dense  about  To  the 
center  of  the  transformation.  is  chosen  to  reside  near  the  peak  of  the 

electric  field  amplitude.  This  transformation  allowed  us  to  obtain  with  60 
computational  time  points  an  accuracy  obtainable  only  with  150  or  better  real 
time  mesh  points.  This  saves  a  significant  amount  of  C-P  time  for  multi¬ 
dimensional  equations.  There  are  penalties  associated  with  the  use  of  this 
procedure.  In  our  case  at  T  *  1  which  corresponds  to  r  near  infinity  the 
inverse  transformation  needed  to  calculate  >T  tends  to  infinity.  Thus  informa¬ 
tion  gathered  near  T  *  1  is  necessarily  regarded  with  skepticism.  In  the 
normal  situation,  the  field  has  decayed  to  zero  at  this  point.  The  number  of 
mesh  points  at  T  *  1  and  T  =  0  is  very  sparse.  At  T  =  0  the  field  is  always 
zero  or  near  zero.  Also,  the  field  at  T*  1  cannot  effect  the  field  at  T  1  -  .*  T 
via  causality.  Pecularities  at  T  =  1  are  discarded. 


The  transverse  boundary  condition  applied  to  the  explicit  scheme 
can  play  a  substantial  role  in  determining  the  propagating  field  when  the  field 
gets  large  at  the  transverse  boundaries,  J.rictly  speaking,  a  linearly  polarized 
electric  field  is  inconsistent  with  the  rotationally  symmetric  vector  wave  equa¬ 
tion  despite  the  fact  nearly  linearly  polarized  fields  do  exist  m  laser  <  s.  illan  rs. 
If  the  walls  are  at  infinity  where  the  field  is  zero,  they  play  no  role  and  ' 1  t 
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boundary  conditions  are  physically  realistic.  In  the  absorbing  type  medium,  this  con¬ 
dition  is  always  satisfied  even  though  the  walls  are  at  n=  1cm.  In  an  amplifying  type 
medium,  the  field  at  the  walls  is  quite  substantial  after  a  few  amplifying  lengths  and 
play  a  large  role  in  determining  the  overall  t -field.  Mathematically  this 
boundary  condition  is  (dpC/C)  =  0.  Numerically  this  is  translated  into 
<*1  =  where  1,  2,  3  correspond  to  the  3  mesh  modes  at  the 

wall  of  the  medium  (o  -  wall  corresponds  to  0j).  This  expression  is  accurate 
to  1st  order  in  Ad. 

The  propagation  of  the  C-field  in  free  space  or  under  the  action 
of  a  focusing  lens  is  handled  by  a  simple  quadrature  integration  in  subroutine 
LENS.  In  free  space,  one  can  obtain  the  solution  of  the  reduced  wave  equation 
by  examining  a  good  table  of  integrals.  The  Green' s  function  propagator  can  some¬ 
times  be  integrated  analytically  with  some  reasonable  input  field.  However,  for 
the  case  where  the  input  field  at  some  z  is  given  only  at  specific  mesh  points  the 
integration  is  easily  handled  numerically.  Care  must  be  taken  to  insure  enough 
mesh  points  are  given  to  handle  large  phase  variations  over  the  aperture. 

By  comparing  results  with  a  known  integrable  case  (Gaussian  input)  we  ex¬ 
perimentally  determined  how  many  2TT-phase  cycles  are  allowed  to  obtain 
reasonable  results.  The  finite  difference  equations  for  the  coupled  equations 
are: 
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l  correspond*  to  longitudinal  position  z  ^  -  t  A  z 
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To  fvirther  understand  the  numerical  procedure,  the  flow  charts  of  the 
main  and  integration  subroutines  are  given  in  Figs.  Bl,  B2,  and  B3 
with  additional  descriptions.  The  main  subroutine  acts  as  an  executive 

supervising  routine  which  controls  the  flow  of  computation.  Its  iwui 
functions  are  to  calculate  the  electric  field  at  z  =  Az  , 

refresh  memory  allocations,  call  input  and  output  functions  and  terminate 
execution.  The  integration  of  9,  W,  and  t  requires  three  space  or  time 
levels.  Thus,  to  obtain  l  at  z  =  2  A  z  will  require  C  at  z  -  0,  and  z  -  Az.  This 
method  cannot  be  applied  in  determining  t  at  z  =  Az.  The  main  subroutine 
calculates  l  at  z  =  Az  using  a  simple  iterated  Euler  finite  difference  approxi¬ 
mation.  This  calculation  of  l  at  z  =  A  z  uses  subroutine  "SUM"  but  the  logic 
for  this  one-time  only  calculation  is  not  included  in  the  flow  chart  for  SUM. 
Subroutines  "LYGRE"  and  "CRAH/M"  perform  almost  all  input  and  output  func¬ 
tions,  respectively. 

The  purpose  of  the  flow  charts  is  to  give  the  reader  the  overall 
flow  of  the  logic.  It  is  not  a  detailed  road  map.  The  detailed  flow  chart  would 
include  various  straight  forward  programming  techniques  to  save  memory 
space  and  time.  The  logic  for  conversion  of  to  gT  fT )S T  of  the  stretching 
computational  space  is  not  given.  Also  there  is  no  need  to  have  9  and  W  as 
two-dimensional  arrays  indicated  in  the  diagrams.  Making  them  one  dimen¬ 
sional  also  requires  the  inclusion  of  a  special  output  format  instruction.  Such 
a  detailed  flow  char:  would  be  unnecessarily  complicated.  /  listing:  <  : 
program  with  copious  comments  is  included  below  for  the  interested  reader. 
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Fig.  B2  -  Integration  Subroutine  SUM 
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Fig.  B3  -  Integration  Subroutine  SUM 


Stability  analysis; 


In  real  space  time  the  reduced  wave  equation  is 

a.e  +£a  e-ic.^e  =  &  .  <b.  3) 

t  n  z  to 

No  matter  what  scheme  is  used  to  numerical  integrate  this  equation,  one  should 
only  gather  information  in  the  light  cone  of  the  input.  There  is  the  obvious 
physical  restriction  that  we  cannot  numerically  estimate  the  electric  field  at 
t  +  At,  z  + Az  from  the  fields  at  t,  z  if  (cAt/Az)  >  1.  If  this  is  violated,  it  means 
the  scheme  is  systematically  introducing  unphysical  data  which  will  lead  to 
meaningless  results  (usually  an  instability).  Thus  a  stability  analysis  of  any 
scheme  applied  to  Eq.  (B.  3)  will  have  the  conditional  requirement  (c'nHAt'Azl  = 
const  <  1  in  the  limit  At,  Az  -•  0.  A  practical  way  to  avoid  this  instability  is  to 
transform  to  a  frame  moving  at  velocity  c/n.  This  procedure  has  been  demon- 

(4 ) 

strated  by  Hopf  and  Scully  in  their  numerical  work.  This  transformation 
has  two  important  benefits  from  a  numerical  point  of  view.  The  transformation 
to  this  co-moving  frame  reduces  by  one  the  number  of  independert  variables, 
thereby  sav'.-g  considerable  calculation  time.  Another  useful  result  of  this 
transformation  is  that  it  automatically  insures  physically  realizable  discretized 
integration  in  the  z,  t  plane.  The  transformation  is  defined  as 


For  any  positive  At  including  the  limit  A"--*  0,  the  ratio  (cAt  'Az)  is  always  greater  than 
or  equal  to  unity  independent  of  the  size  of  An  =  Az.  Thus  the  stability  with  respect  to  **.- 
is  built  into  the  integration  scheme.  A  drawback  of  this  type  of  transformation 
will  occur  when  trying  to  handle  pulses  with  velocities  substantially  different 
from  c/n.  Since  only  a  finite  number  of  time  mesh  points  will  be  used  to 
decretize  the  time  domain,  it  follows  the  time  window  in  the  co-moving  trame 
is  finitely  restricted.  If  a  pulse  has  a  velocity  less  than  c/n,  after  passing 
through  a  calculable  distance  z  into  the  media,  the  pulse  will  eventually  drift 
..’it  of  this  fixed  time  window.  This  could  be  compensated  for  by  using  a 
variable  velocity  co-moving  frame.  However,  the  code  becomes  more  complex 
and  time  consuming.  Fortunately,  for  the  transient  type  problems  \\t  have 
analyzed,  this  difficulty  does  not  restrict  the  numerical  procedure.  In  appr.  ach¬ 
ing  steady  state,  the  problem  becomes  more  serious. 
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The  numerical  scheme  that  wt  have  used  effectively  treats  the 
coupled  equations  like  two  coupled  ordinary  differential  equations.  ’ATnle  r 
is  held  fixed,  the  field  equation  is  marched  in  z  by  a  Az,  and  when  z  is  held 
fixed,  i3  and  W  are  stepped  in  time  T.  The  stability  analysis  is  applied  to 
each  of  these  processes  separately  and  the  combined  criterion  is  applied  to 
the  coupled  equations.  It  can  be  shown  that  for  these  equations,  this  result 
is  valid  as  well  as  easy  to  apply.  This  result  is  more  apparent  once  the 
analysis  is  demonstrated  for  the  two  separated  calculations.  A  complete 
description  of  the  methods  of  numerical  stability  is  given  in  Richtmyer  and 
Morton.  The  stability  check  is  applied  locally  to  the  linearized  equation. 
One  of  he  separated  integrations  carried  out  is  for  the  case  g  =  0.  The  re¬ 
duced  wave  equation  becomes  a  hyperbolic  diffusion  type  equation  with  an 
imaginary  diffusion  constant  "iCj.  "  After  a  Fourier  series  expansion  of  the 
discretized  dependent  variables,  the  finite  difference  equations  can  be 
assembled  in  the  form 


Vn  +  1  r  A  Vn 


rn+l  . 


where  V  is  the  vector  representation  of  the  dependent  variables  at  the  mesh 
point  zn+i  and  A  is  the  so-called  amplification  mat  rix.  The  VonNeuman  neccss it\ 
condition  for  stability  requires  the  amplification  matrix  to  be  bounded.  This  generally  re¬ 
quires  its  eigenvalues  to  have  absolute  values  less  than  or  equal  to  unity.  The 
detailed  calculations  and  results  are  shown  in  the  next  section.  The  other  effec¬ 
tive  ordinary  differential  equation  will  be  the  case  where  Cj  0  (the  one-dimen¬ 
sional  problem).  The  results  of  the  3-level  predictor-corrector  scheme  (levels 
n  -  An,  n,  n  +  An)  requires 


2Anc 


~2  -  °-4q 


1  in  our  calc,  we  used  v  -  0.  33 


and  no  condition  on  A“  ,  A-  other  than  they  be  consistently  small.  Via  numerical 

experimentation,  we  have  found  the  cut-off  value  of  \  to  be  •  "0.4  Thu  a  tlu 

c 

local  stability  check  comes  reasonably  close  to  the  actual  value  (as  is  usual!, 
the  case).  Global  instability  almost  always  results  from  local  instability.  In¬ 
cluded  in  the  next  section  is  a  figure  plotting  the  absolute  values  of  the  ampli¬ 
fication  matrix's  two  eigenvalues  >  i  and  >  for  various  values  of  the  critical 
parameter  '<  (2Anc  ]  /  (Ao  )“.  Again,  stabilitv  re  quires  h  |  <  1. 
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Stability  check  for  p-i  integration  scheme 


-■ase  ia-  g  -  0  Reduced  wave  equation: 

The  3-level  finite  deference  solution  for  l  is: 

'3"‘‘ =  +  *"♦'*  * 5? te^ +  'r  “jVi  *  ~z*U  *  e"-2'jntI 


where  n  -  T]  and  j-y.  The  Fourier  series  substitution  is  simply  accomplished 
by  letting 

X£e±iHAy  replace 

i  ±  1 


l  eii5tAV 


for  g"  , 

j±  1  ' 


sn  =  e11 "  1 


The  resulting  amplification  matrix  with  this  substitution  is  (2x2) 


k  where 

-1  X 


A  U 


3  r  X  "  \  +  'l  f Y  sin2  (~L)  +  2y 2  sin4 


(^f  ) 


b  =  ■  I  +  i  2  sin2(JTZ)  • 

The  characteristic  equation  is  a>  +  b  -  0.  After 
the  absolute  value  of  its  roots  are: 


some  straightforward  algebra 


\c  -  i  f  y[  +  \/x2  +  y2 


2x,(^  +1  /x24v2)1/2  t2v *  *  I  /  2  i  ..2  1/2 


T2  r  2 


2ylf“  i  +  i  S  i  y  )  W“  t 


where 
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1  .  2  .  4,kAv, 

ij  =  2  +  2y  sin  ' 


Vl 


5  .  2 ,  v  Ay  v 

=  -  -  y  sin  (— p) 


X  =2-fv2Sin4(^)  +  4v4Sln8(^f» 


and 


y  =  -  |  v  sin2(il|^)  +  10v3  sin6 


2  T  '  2  '  ‘  ”  ’  .  2 

For  all  real  values  of  y  and  H  .  the  sin  is  restricted  4o  the  range 


0  <  sin2(  1 

The  absolute  values  of  the  roots  |  X±  |  are  plotted  over  this  range  in  Fig. 
B.  4. 

Case  2.  c  =  0  coupled  one-dimensional  linearized  system  equations  arc 


=  -ce+gP 

3  P  =  -a P  -  a£  +  bW 

T 

aw  =  k  -  flW  -  c6  -  dP  . 

T  o 


iiripr  writing  the  complete  three-step  complete 

.  i  t  x t  nl  .  ,  .  ±in A* 

+  i  w  +  and  P  ,  and  substituting  hxe 

n  1 1  ’  n+1  n+1 


discretized  equations  for 

.  r  4 1  .  ,  . 

lor  ixn_.  j  yield, 


).W[/“At  -  |  -f”  |  4  4t[-  •  •  I  -  0 

.  ,  -in  At 


and 


>  s  =  e 

Where  L  corresponds  to  V  and  n  corresponds  to  t.  Letting  At,  Ar  -  f 
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±  v 
2 


<  1 


x±  =  m 


TO)f2 


This  verifies  the  stability  of  r,  n  integrations.  We  note  that  no  mutual 
strictions  are  placed  on  At,  An  as  noted  in  text  on  physical  grounds. 
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